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ABSTRACT 


The Navy Remote Ocean Sensing System ( N-ROSS ) satellite is being developed 
to supply accurate data on ocean parameters for fleet operations. A Low Frequency 
Microwave Radiometer ( LFMR ), a large flexible reflector attached to an angled 
flexible boom, is a sea surface temperature sensor on this satellite which rotates at 15 
RPM. The dynamic interaction between the reflector and the boom, and the effects of 
the reflector orientation and flexibility on the pointing error of the LFMR during a 
Spin-up procedure are investigated by performing dynamic simulations. Dynamical 
equations of this flexible multibody system are formulated using Kane’s method. 
Efficient computer simulations were achieved by developing a FORTRAN program 


and using Dynamic Simulation Language ( DSL ). 
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I. INTRODUCTI 


A. BACKGROUND 

The Navy has great interest in obtaining accurate, real-time data on various ocean 
parameters. Presently, research is being conducted on the merits of using a satellite system 
to collect this data. The Navy Remote Ocean Sensing System (N—ROSS) satellite has been 
developed to achieve this goal [Ref. 1]. This satellite scans the earth’s surface and will 
provide the Navy with such parameters as wind speed, wind direction, ocean temperature, 
ice edge detection, wave height, and ocean photography [Ref. 2]. 

One of the sensor systems to accomplish this goal is the Low Frequency Microwave 
Radiometer (LFMR). The LFMR consists of a large reflector dish connected to an angled 
boom. The LFMR scans the ocean surfaces and records the temperature. Increased scan- 
ning area is accomplished by spinning the reflector at 15 RPM. A rigid electronics box is 
attached at the base of the LFMR boom. Therefore, the rotating LFMR system consists of 
the flexible boom and reflector and the ngid electronics box. Considering the nature of the 
system, accurate analysis and suppression of the vibration and deflection of the LFMR is 
essential for the satisfaction of the pointing error requirement of the system and for the 


accomplishment of the mission of the N-ROSS satellite[Ref. 2]. 


B. PROBLEM STATEMENT 
A three-dimensional dynamic analysis has been performed using a model that consists 
of a flexible angled boom with a point mass at the tip [Ref. 3]. Lagrangian equations of 


motion were employed to achieve a good dynamic simulation. 
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In this study, the LFMR is analyzed as a multibody system. A multibody approach 
will provide a means to analyze the dynamic interaction between substructures of a compli- 
cated structure with greater ease. In dynamic analysis of a multibody system, maintaining 
computational efficiency is a key factor and “The emphasis of researchers working with 
multibody systems has been the expanded generality of mathematical models and the for- 
mulation of equations of motion that are amenable to computer solution.” [Ref. 4] 

Therefore, it is the intent of this thesis to apply an efficient multibody dynamic analy- 
sis to the LFMR system. The LFMR is broken down into two bodies; the first body is the 
flexible boom and the second body is the flexible reflector. Kane’s method is used to for- 
mulate these equations of motion, “since it combines the computational advantages of both 


Newton's laws and the Lagrangian formulation.” [Ref. 4] 


C. THESIS OUTLINE 

In Chapter II, the analytic model of the LFMR used for this thesis is described. The 
dynamical equations are formulated using Kane’s method. This method is briefly 
explained prior to the formulation. 

Chapter II contains an explanation of the three computer programs used to solve the 
equations of motion. NASA Structural Analysis (NASTRAN) calculates the mode shapes 
and natural frequencies for both bodies. A FORTRAN program reads the data output from 
NASTRAN and calculates the time constants. The time constants are applied in a Dynamic 
Simulation Language (DSL) program to solve the simultaneous differential equations. 

In Chapter IV, the simulation results are presented to investigate deflections and rota- 
tions at various points on the LFMR. Three parameters are varied: (1) the flexibility of the 
reflector; (2) the in-plane orientation of the reflector; and (3) the out-of-plane orientation of 


the reflector. The effects of these changes are compared. 


J.) 


Chapter V presents conclusions and recommendations for application and extension 


of this work. 


WZ 


MULATI F ATION 


A. MODEL DESCRIPTION 

The configuration of the LFMR is shown in Figure 2.1. The LFMR consists of a 
boom mounted on a ngid electronics box. A complex reflector is attached at the tip of the 
boom. The boom hinge is to be rigid after the deployment of the LFMR, allowing no rela- 
tive motion between the upper and lower booms. The LFMR system rotates at 15 RPM 
about a spin axis fixed in a local reference frame of the N-ROSS satellite. The spin axis 
always points at the earth’s center while the satellite is in its orbit. Therefore, the gravita- 
tional force and the centrifugal force are assumed in equilibrium. This observation allows 
the spacecraft’s local reference frame to be considered the inertial (Newtonian) reference 
frame and the LFMR to be in a zero-gravity environment. 

The model of the LFMR used in this thesis is shown in Figure 2.2. The angled boom 
was analyzed as one body and the reflector as the second body. Body A represents the 
electronics box, body B is the angled boom, and body C is the reflector dish. 

The @;,a2,a3 coordinate system rotates in the inertial reference frame, with az 
coinciding with the spin axis. The C;,C2,Cz coordinate system is the local reference frame 
of the reflector and its origin is at the connection point between the deflected boom and 


reflector (point h). 


B. EQUATION FORMULATION 


1. Kane’s Dynamical Equations 
The equations of motion for the LFMR were formulated using Kane’s method 
(Ref. 5]. This method is a Lagrange form of D’Alambert’s principle, which is equivalent 


to Newton’s law cast into a different form. Newton’s law for a differential element is 
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(Eqn. 2.1) df - adm =0 
where 
df - differential force 
a - acceleration of differential element on body K 
dm - differential mass 
Kane’s method introduces the generalized active and inertia forces. The 
dynamical equations are 
(Eqn. 2.2) F.+Fo=0 r=t,..., total degree of freedom 
F . represents the generalized active forces and ™ is the generalized inertia 
forces. Fis the sum of the dot product between the partial velocity and the differential 
force as shown in equation 2.3. 
(Eqne2s3) = Jy. suf, 
= is the sum of the dot product between the partial velocity and the differential 
inertia force as shown in equation 2.4. 
(Eqn. 2.4) Jy. Wee cio 


The partial velocity V, comes from the definition of velocity where 


V 
(Eqn. 2.5) v= Dy. q r=l,..., total degrees of freedom 
r= | 


q, is the time derivative of the generalized coordinates. 
For this model, there are three types of dynamical equations. Equation 2.6 is 
derived from the large rotational motion of the system. Equation 2.7 is derived from the 


small boom motion. The number of equations of this type depends on the number of 


modes (V ) used to characterize the boom motion. Equation 2.8 is derived from the small 
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reflector motion. The number of equations of this type depends on the number of modes 


(V>) used to characterize the reflector motion. 


(Eqn. 2.6) Ei Ey =e0 
eon: 2) F yj Bs F 45} = 0 Jae 
(Eqn. 2.8) F 1+vy+k F Bente = 0 Kole. V2 


2. Generalized Coordinates and Generalized Speed 
The first generalized speed U , is defined as follows: 
edi 2.9) UE ee ae a el 
where at is the angular velocity of body A in the inertial reference frame. Subsequent 
generalized speeds are defined from the generalized coordinates q,, qj which are model 
coordinates of bodies B and C, respectively, as follows: 
(Eqn. 2.10) ine =q: (ae) 
Tey) = ol (j = 1, ..., Vo) 
where V; and V> represent the number of modes used to describe the displacement of body 
B and C, respectively. Thus, there are 1 + V; + Vo generalized speeds. 
3. Angular Velocities 
The angular velocity of body C with respect to body A is the time derivative of 


the small rotation of the hinge point H. Point H is the connection between body B and 


body C. This angular velocity can be represented using matrices as follows: 


foe - | ¥3-Y2 Q -Y3 2 0 -P3 Y2 
(Eqn. 2.11) “w =lc} [c] =| -¥3 | Y3 O -Yi f=} ¥3 0 -¥ 
Yo-Y, | muy ye “YoY, 0 


V4 
, 
where ¥; = > Y. (H) Uys, and ¥, is the small angle displacement for the i! 


degree of freedom and the k'® mode at point H. Therefore, the angular velocity is 


Is 


AC, 
(eGinee ei) (iC 1 tio Comme aes 


4. Velocity and Accelerations 
a. Velocity and Acceleration of Body A 


From the generalized speed definition, the angular velocity of Body A is 


NA _ 
(Egne 2:13) ee = Ue 
The angular acceleration is merely the time derivative of the angular 


velocity, because the az axis is fixed in the inertial reference frame. 


NA. 
(Eqn. 2.14) (O66 Sau Pek 

The velocity of A is defined relative to point Q, the origin of the a;,a9,az 
coordinate system. Point Q is fixed in the inertial reference frame. Since body A is a rigid 


body, the velocity is merely the rotational velocity of the mass center A*. 


ea NA 
(EGiie 2s ie) Vv =(q@ x -ba)) = -DU sae 


where D is the distance from point Q to the mass center A*. 
Similarly, the acceleration of body A is the rotational acceleration com- 


prised of tangential and normal components. 


NAX NA UNA NA 
(Eqn. 2.16) a = WW x(W * -DapiceeaGra, 


2 
Dugmapeeeel, ao 


The 2;,a2,a3 coordinate system was chosen so that the mass center of 
body A moves in plane motion, simplifying the equations. 
b. Velocity and Acceleration of Body B 
The position of an arbitrary point P on body B is described in the 


a;,a>,a3 coordinate system. 
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V4 
Be x 
(Ben. 2.17) im Secunda i. az + >>; 4 ay 
i= | 


V1 y V1 ; 
: 2 qi 42° a qj 23 
j= = 


The position vector a is the sum of the undeformed position vector and 
the summations of the translational mode shape deformations. P,,P>, and Pz describe 
the undeformed position. 

The position of point P changes with time, therefore the velocity of point 


P is the time derivative of the position added to the rotational velocity. 


NP B.P NA BP 
V R 


(Eqn. 2.18) *+q@ x RB 


V4 ‘ V1 y Vy : 
: a; Ujey a1 * 2° Uy+; d2 * 2X; U je; as 
ls |= | 


|= | 


+ 


val V4 
Uy Pit Doig, de> Uy Pi+ DD o%q, gy 
)= | l=] 


The acceleration of point P is derived from differentiating the velocity 


with respect to time. 


(eqng 2.19) x = 


ky 


V 
2U | aor bP] > a; * 


|=] 


Th ai 54% qj uf ie Sails 


}=] 1= | 


V1 
Zl) To une +o U4; fol, 
|= | 


l=] > 


V1 


i > $5 U i+) Bs 


|= | 


c. Velocity and Accelerationof Body C 
The position vector of an arbitrary point X on body C is described in the 


local coordinate system C, analogous to the representation of a point on body B. 


(Bene. 20) IB, = “ae > 1G, Si” | gee * > 


The velocity of point X in the inertial reference frame is the sum of the 


velocity at the hinge point H in the inertial reference frame and the rotational velocity with 


the time derivative of the position calculated in the local coordinate system C. 
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(Eqn. 2.21) Vs eerie) + R 
where 
NC 
W@W =U; ast ¥, C1 + Yo So * V3 G3 
NH ai x ua 
v = > (H)Uy. 47 U | Boe: b* (H)q, fol = 


Ot ea 


b5 (Hu. tU; Bai d $;(H)q, alp © 
i= 1 


V1 
2X PiCH)uj.1 83 
|= | 

Thus the velocity of point x is 


V1 V 1 
NX x 
(Eqn. 2.22) 'V° = 2 $i Huier ui] Hot Oy CHG; | Par 
= 


f=] 


Loud ao + 


i 
DX oi H)ujye tua] # 
=| 


iT noe 


V1 V> : 
|= ] j=l i= 


Me 


a 2 3 5 v2 xX 
| eT 2° Do) Y 2 xt 2 $54 Ca 
Nie l= j= 
Vo _ 
nou 2 915 ds. 
le 


A third nonorthogonal coordinate system (d) is defined for computing the velocity 
and acceleration of body C more easily. 


ES = 
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The following tables depict the relationship between the three coordinate systems. 


ey 
d> 
Q3 


gy 
Q> 
G3 


gy 
Q> 
Gs 


(not orthogonal) 





C1 
C2 
C3 





The acceleration of point X with respect to the inertial reference frame is 


expressed as the hinge point H acceleration and the acceleration relative to the hinge point 
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NX 
(Eqn. 2.23) a = @ + ox < wBh 4 G0 Gee eee) 





where 
AC 
NC NA “d(w) NA ac 
ORS OS dt pO 
=U, as + Y, oS) + Vo Sot Wz cz + Uy ¥) Gy 
=U, Vo do + Uy ¥3 ds 
and 


Vv 4 V1 

NH ; ZL 

a =4-U, Har Doi Wa us Hie Leia 
i= i= 


Vv 4 val 
~2Uy 2 our 2 ids. a) * 
i= i= 
V1 . , V1 y 
U | Hit 2D >; (Hq, el Hot >, O5(H)q, 
l= | l= 1 


V1 V1 
-2uy diye DV oy(Wuy.; P ae * 
i= | |=] 


V1 : , 
l= | 


Te 


The acceleration equation for body C was too unwieldy for higher terms 


to be considered. Therefore, only first-order terms were used to produce 


We 
NX NA 2 Z a Tes 
omee 24) a= a | ou ey pu Ao tu x3+ 2 O54, a3 
Ie) 


V2 
, | 2 cee 
+] Yoxs-VsxotuyPyxg-Uy] Xt DO jG) P¥suixy + 
i 


V 


N 
<< 


2 
x. | , | 2 =U. 
® Uyey, <7 | Sit | —¥i%3t Vax, t Ui V2x3-U 4! Rot d o79, 
! j=] 


w@ 
ql 2 7 eS 
— ¥3UyX2+ DD, Sauces Seale cm oe i| Xs > qj 
|= 


——a 
il 
—— 


We 
, Teac ‘ Tae ee 
ss nae oe |) al 19, 
)=I )= 


V2 
’ ; ee 8 
* Uy Yo X3 — Uy V3 Xo + Uy 22 Ytovyep ae 
7 


Ze 


Oi 3%, 


- uy, Yy %3 + Uy Pz Xy d+ 
j= 


<< 
N 
— 
N 


=>, -7 
2Uy (eye lee Cen > | 


Nat © 
i 
—, 


Me 
WS) ee Ue Se See 2 uF) iene nos 


Partial Velocities 


The generalized partial velocities and angular velocities are assembled in Table 


2.1 and Table 2.2. These values and the acceleration equations are the basis of the formu- 
lation for the generalized inertia forces. 


TABLE 2.1 
PARTIAL VELOCITIES OF BODY A AND BODY B 


r Ww’ Ve ye 
V1 . 
“Dae as os gi | 1 * 2919 oe 
|= 
Z 
a) + $5 a2 * > | as 
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ABE ES 2.2 
PARTIAL VELOCITIES OF BODY C 


e S 
r w, ver 
V1 y Vy g 
23 1) oa > $5 (H)q, ajt | Hyt » >; (H)q; le 
l= 1 |= | 


y (Hor + |b) (H) ay + 4 CH) ao + $5 (H) as + 


2 2 <2 Zz 3 < ry 
Yj CH) cot HY; st 2s Pia =e ot de Pid ©] 


we ae | “2 Tes 
ae | it 2s 4a a 3" 2s Pak <2 
M2 2 
Te 2 re 
Yj a" 2s Pia a Es ex 
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6. Generalized Inertia Forces 
a. Generalized Inertia Force of Body A 
The generalized inertia force for body A is comprised of two parts. the 
first is due to the moment of inertia of body A’s mass center about the az axis. The second 
is due to the inertial force of body A’s mass center. 
(EqneZ 25) AR = w (~J. ai — Ma, a m 
where Jz is the moment of inertia of A* and Ma is the mass of body A. 


A A* 
Using the values of @), and V , from Table 2.1, the generalized inertia 


force of body A is 


(Eqn. 2.26) AF, = J3 Uy - Mm, D2 u, 


There is only one equation since body A is rigid and displays only rota- 
tional motion. 
b. Generalized Inertia Forces of Body B 
Body B is a continuous elastic beam with a mass per unit length of p. 


Thus, the generalized inertia equations are an integration along the length of the beam 
x NP 
(EqueZz 7) 2a =- f (ue -a) pdx 
0 


The result is an equation for the large motion and one equation for each 
deformable degree of freedom of body B. 
The large motion equation incorporates only first-order terms, while the 


small motion equation neglects third-order or higher terms. 
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(Bans 2°25) om = f Uj pi +P5+2P, x4 Sue ey, 5 4! qj 
Q |= | 


V4 


Es | Bi {Pitt ati ZU Buf $5 +P od 


< 


Ddx 


ae 
—S$_— 


AF of 
|=] 


il $y eet -$7| P y+ SN 
Q |= | 


Er sil ee aac 


iei-atet p20, Zora {tte-bht foo 


c. Generalized Inertia Forces of Body C 


—m « 


Body C was modelled with lumped masses at 63 grid points. Therefore, 


the generalized inertia equations consist of a summation of the acceleration dotted with the 


generalized partial velocity times the mass 


Ly 
is NX 
(Eqn. 2.29) io »y (si a J 


Pa) 


There are | + V; + Vo equations of motion for body C. Again, the 
large motion equation has only first-order terms and the subsequent small motion equations 


neglect third-order and higher terms. 


(Eqn. 2.30) 
a 2s Ze st a <¢. 
y= Duy FH ytHoexytxgt2CHi ex) Dia; | W$ja + 
v=G j= | j=l 


V4 One 
2 (He + xo) | LPi(Hay+ DFG [+ 2xiHy + 2xoHe - 
Bi ie 


ZY 1X3 (H> a Xo) at 2YoX2 (H a on) ar 2Y3 (Hox a H Xo) a 
va x y 

2D, Uae |-#000Hee nad FONE “Ph rstHa xD = 

j=] 


"(Hot xo)*P2[x (xy +H Deraberhel ae 


<2 7% re! V4 
2 Yisivs | Fitertad Fit rmd 2uy Db 
iE z 
Tarn @HN Hare + X3 y, (Ho r X95) _ 


x3 yi (Hy * X4) yy China 1) : 
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Bry [4 ne $; ver] fs 


. L2 ar res, 
Saas ath > (H) Y1X3-%3X1-H2-X2- DP} (H) qj - IK 
|= | K= | 


x=0 


V Zee 
p) (H) | Yoxs YswarHirxrs Wei H)a) + 2 64a 
|= = 
V1 
1G Xo | - = Silly Me ings ales teea Dia + D Fea 


Vv | ae 
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(Eqn. 2.31) 
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7. Generalized Active Forces 
The generalized active forces are derived from the forces exerted on the differ- 
ential elements. Only the internal elastic forces and the external forces were considered in 
this formulation. Rotation of the individual beam elements was neglected. 


The equation for the generalized internal forces is of the form 


Lo 
(Eqn. 2.32) F = j y,8- dfB+ div. afc 
6 x=0 


where df consists of internal axial and shear forces. 


The shear force portion of df for the beam is 


ae acy 
(eda. 2.33) Ol shear = By2 | Elaya | dx ae 
V1 , 
where Y = Xd: qj 
ia 


The axial force portion of df for the beam is 


ae 
dtayia, = AE Axe dX ay 


V4 
x 
where U = do; qj: 
J=I 
When these differential forces are dotted with the generalized partial velocity 


and then integrated by parts, the remaining terms are simply the forces due to strain energy. 


The generalized internal forces are reduced to 
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| zB 
(Eqn. 2.34) F eee OU 


The external force on the boom was the torque applied about az axis at point Q, 
the fixed end of the boom. This torque was dotted with the generalied partial angular 


velocity to produce 


E 


For the reflector dish, it was assumed that the strain energy would account for 
the internal forces and no external forces were applied to the reflector. 


Thus, the generalized active forces for the reflector are 


| 0 oe 
(Egn. 2.36) F jay yk ——a QD j My a, 


8. Final Equations 


The equations of Sections B.6 and B.7 were gathered to obtain the following 


equations: 

(Eqn. 2.37) Afi + Bro + Cha + FE = 0 

(Eqn. 2.38) BF ey + BF ey + CF yj =O j= ly 
(Eqn. 2.39) CF avis + OF uae =0 k = 1,...,V9 


The orthogonality condition for the boom is: 
f x |x y iy ae 
(Eqn. 2.40) ; Pi >] * PEP] t Pi Pj POX * 


Lo ’ 
> sieogeoo-Blenpo diane |e 


x=0 
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asi) 
= 0 = j 
and for the reflector: 
=2 Tor een ere 7 2 =, 7 
x= 
=a | = | 


These conditions were incorporated into the final equations to reduce the 
number of terms in the final equations. 
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I. COMPUTER IMPLEMENTATION 


A. INTRODUCTION 

Computer implementation to solve these equations of motion involved three stages. 
Modal analysis was conducted using the NASA Structural Analysis (NASTRAN). 
NASTRAN is a general-purpose digital computer program for the analysis of various 
structures [Ref. 7]. NASTRAN provided values for the natural frequencies, mode shapes, 
and generalized modal masses for the two flexible bodies. 

The data was transferred and read into an IBM/VS FORTRAN program. Constants 
were calculated in this program, which then were inserted into the initial stage of a Dynamic 
Simulation Language (DSL) program. DSL is a digital simulation for continuous systems 


[Ref. 8]. DSL was used to solve the simultaneous differential equations. 


B. MODAL ANALYSIS 

A multibody dynamic analysis requires special consideration of the boundary condi- 
tions. For the boom, one end was fixed and the other end was free with a tip mass 
attached. The tip mass was equal to the total mass of the reflector. The mode shapes, 
natural frequencies, and generalized masses were calculated using 17 grid points (Figure 
A.6). Figures A.7, A.8, A.9, and A.10 depict the first four modes of the boom. 

The reflector was represented by 63 grid points (Figure A.1). Various links between 
grid points were constrained from rotating and the hinge point (Grid 63) was fixed. 
Figures A.2, A.3, A.4, and A.5 depict the first four modes of the reflection. 

Ten modes were generated by NASTRAN for each body. A preliminary analysis 


between two, three, four, and five modes for each body was conducted. Four modes was 
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chosen for the final analysis because minimal changes occurred between the four- and five- 
mode cases. 

Modified Given’s method (MGIV) was used to conduct the modal analysis of the 
boom and the reflector. The generalized mass was normalized to equal one simplifying the 
equations. 

The natural frequencies for the boom and two cases of the reflector are tabulated in 
Table 3.1. In the results section, the effects of reflector flexibility are discussed. The 
NASTRAN programs are shown in Appendix E. 

TABLE 3.1 
REAL EIGENVALUES 


Boom Reflector 1 Reflector 2 
Natural Natural Natural 
Mode Frequency Toles). Frequency 





(cycles) (cycles) (cycles) 


C. EVALUATION OF TIME CONSTANTS 

Two files were transferred to the IBM/VS computer. Both files contained ten modes 
of data for one of the bodies. All ten modes are read into a FORTRAN program, but the 
number of modes used to calculate constants was varied by a parameter statement in the 


beginning of the program. This strategy provided easier formatting to read in the data, 
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while still generalizing the program and making the computations more efficient. The pro- 
gram is shown in Appendix E. 

The equations of motion were systematically evaluated to identify terms that depended 
on position alone and were repeated several times. These terms were isolated from the 
time-dependent terms. Subroutines were written to calculate these position terms and were 
called immediately after reading in the mode shape data. The main program consisted of 
multiplications and additions between system constants and the subroutine outputs, since 
the subroutines carried out the required position integrations or summations. 

The reason for having separate programs is that all the position calculations can be 
completed beforehand. Thus, the DSL program need only be concerned with time step 
calculations. The number of constants and their dimensions depended on the DSL imple- 


mentation and will be discussed in that section. 


D. DSL IMPLEMENTATION 
With the separate FORTRAN program to evaluate the time constants, the DSL pro- 
gram is relatively simple. The equations of motion are put in the form 
[Alu=8 
where [A\] is a matrix of the coefficients of the acceleration terms (the left-hand side of the 
equations of motion), while B is a vector of all other terms (the mght-hand side of the 


equations of motion). 

These coefficients were supplied by the FORTRAN program. The A matrix and B 
vector were assembled in the derivative section of the DSL program. Two subroutines 
from LINPACK were called to decompose the A matrix and solve the matrix equation. 


DGEFA uses gaussian elimination to decompose A and DGESL solves the matrix 


equation. 
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DSL provided many alternatives to solve a time-stepping problem for simultaneous 
differential equations. A Runge-Kutta fifth-order method with variable step size was 
selected because it was self-starting, stable, and accurate. Although other methods may 
have been more efficient computationally, accuracy was a primary selection criterion. The 
program is shown in Appendix E. 

The results were printed out in a file and plotted using the TEK618. Plotting was 
accomplished via the GRAFAEL command, which 1s a particularly powerful accessory of 


the DSL program. 


aN, 


IV. RESULT 


Three parameters of the LFMR system were varied in the computer simulation to 
investigate the effects of these parameters on pointing error of the LFMR system during a 
spin-up procedure. Deflections and slope changes at the boom tip (grid point 63) and grid 
points 1 and 4 of the reflector were compared for this purpose. The locations of these grid 
points are shown in Figure 4.1. The first analysis compared the effects of changing the 
reflector flexibility, while the boom properties remained the same. The second analysis 
compared vibrational amplitudes of three cases with different in-plane orientation of the 
reflector. These orientations were accomplished by lining up the a;,a2,az3 coordinate 
system with the C;,C2,Cz coordinate system and rotating the C; and Cz axes about the a> 
axis to the desired angle. The third analysis compared the case of a particular in-plane 
orientation (-155° case) with the case of the same in-plane orientation plus slight out-of- 
plane tilt (+5°). The out-of-plane orientation was achieved by rotating the reflector about 
the az axis. 

The same spin-up procedure was used in all computer simulation runs and the torque 


applied to the LFMR system during the spin-up procedure is shown in Figure 4.2. 


A. MATERIAL PROPERTIES AND PARAMETERS 


The material used to model the boom and reflector is Isotropic Graphite-Epoxy 


Composite (T300/5208 (0/90/45/-45).) [Ref. 9]. The inertia property of the electronics box 
was obtained, assuming it to be a uniform two-foot cubic body. The properties of the 
boom and reflectors are presented in Table 4.1. Reflector 1 is the baseline design of the 


reflector used throughout the analyses. Reflector 2 is a more flexible model of the reflector 
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and was used in the reflector flexibility analysis. The geometric parameters for the boom 


and two models of the reflectors are given in Table 4.2. 


TABLE 4.1 
BOOM AND REFLECTOR PROPERTIES 


Modulus of elasticity 10.1E + 6 psi 10.1E + 6 psi 10.1E + 6 psi 
(E) 


Poisson’s Ratio 
(v) 


x-Sectional area 1.1868 in? 1425 in .1425 in 
(A) 

Inner diameter 

Rod element 2.5 in 1.04859 in .8177 in 
(d,) 

Outer diameter 

rod element 3.0 in 1 in 75 in 
(d,) 

Thickness of rod element .04859 in .0677 in 

Mass per unit length 1.522 E-4 1.522 E-4 1.522 E-4 slugs/in? 
(p) slugs/in3 slugs/in3 


Area moment 


of inertia 1.2231 in4 0165 in* 00853 in? 
(A) 


iz) 
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TABLE 4.2 
GEOMETRIC PARAMETERS 


Pe a 


|,: length of lower boom 162 in 
I>: length of upper boom 126 in 
c: angle between |; and az axis 34,5° 
B: angle between |, and lo 98.76° 
M,: weight of electronics box 50 Ib 
Mc: weight of reflector 69/5) |l0) 


B. EFFECTS OF THE REFLECTOR FLEXIBILITY 

Two cases of reflector flexibility were chosen for comparison. As shown in Table 
4.1, Reflector 1 was stiffer than Reflector 2. The difference in the angular velocity is 
shown in Figures B.1 and B.2. Greater amplitude oscillations occur after the applied 
torque is removed for the more flexible reflector. This effect is due to larger vibrations of 
the more flexible reflector (Reflector 2). 

The fourth generalized coordinate (Q4(t)) of the boom was chosen to represent the 
reflector flexibility effects on the boom’s generalized coordinates. Figures B.3 and B.4 
depict this boom generalized coordinate for the two cases. They show a definite dynamic 


interaction between the boom and reflector. The more flexible reflector displayed greater 
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fluctuations in the boom generalized coordinate, which translated into greater vibrations at 
the boom tip. The vertical deflections of the boom tip are shown in Figures B.7 and B.8. 

The third reflector generalized coordinate (Q3(t)) was chosen to represent the effects 
of reflector flexibility on the reflector’s generalized coordinates. Comparison of Figures 
B.5 and B.6 shows a negative shift in the graph of the more flexible reflector. The ampli- 
tude of the reflector generalized coordinate is unchanged. 

The rotational vibrations about the vertical axis are compared in Figures B.9 and 
B.10. The more flexible reflector caused greater amplitude rotational vibration at the boom 


tip. 


C. EFFECTS OF IN-PLANE ORIENTATION CHANGES 

As mentioned previously, cases of different in-plane orientations between the reflec- 
tor and the boom were analyzed: Three cases of angle rotations about the a> axis were 
compared: -135", -145°, and -155°. The angular displacement for these three orientations 
is Shown in Figures C.1-3. The angular displacements oscillated after the applied torque 
was set to zero. An increasing amplitude trend in angular displacements was apparent as 
the orientation was a less negative angle. 

The horizontal deflection at the boom tip represents the effects of reflector orientation 
on the boom. Figures C.4-6 show that greater amplitude vibration occurred for the case of 
135° reflector orientation in the horizontal direction. 

The rotations at the boom tip about the horizontal axis are shown in Figures C.7-9. 
Greater fluctuations occurred for the —135° orientation. The larger amplitude results for the 
case of -135° orientation is due to the location of the center of mass of the reflector. The 
radius of rotation is larger for the -135° orientation. Thus, the moment of inertia is 


increased. 
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The reflector orientation also affects the reflector vibration. Figures 3.10-12 show 
the horizontal displacement of grid point 1. Higher frequency components appear in the 


~—155° case, but the greater amplitude vibrations occur in the -135° case. 


D. EFFECT OF OUT-OF-PLANE ORIENTATION OF REFLECTOR 

The case of -155° in-plane orientation of the reflector 1s compared to the case with the 
same 155° in-plane rotation plus a 5° tilt out of plane. The first boom generalized coordi- 
nate (Q ,(t)) was the only boom coordinate affected by the tilt out of plane. The tilt pro- 
duced a positive shift in (Q ,(t)) value as shown in Figures D.1 and D.2. The oscillations 
were same in amplitude. 

The fourth reflector generalized coordinate (Q4(t)) was chosen to represent the 
effect of the out-of-plane orientation on the reflector generalized coordinates. Figures D.3 
and D.4 show that vibrations with larger amplitude appeared for the out-of-plane 
orientation. 

The out-of-plane orientation deflections for the boom and reflector were the most 
affected by the out-of-plane tilt of the reflector. Figures D.5 and D.6 show the effects at 
the boom tip. The first boom generalized coordinate was the predominant component 
resulting in a positive shift in the out-of-plane displacement of the boom. Figures D.7 and 
D.8 show the effects at grid point 4 with the first boom generalized coordinate producing a 
positive shift. The reflector generalized coordinates produced a larger amplitude vibration 


for the out-of-plane rotation. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The purpose of this research was to develop a realistic model of the LFMR sensor 
system. This goal was successfully met by applying Kane’s method to a multibody model. 
Efficient computer simulation was achieved. The programs were general enough to use up 
to ten modes for each body, but the results indicate that superimposition of four modes for 
each body was sufficient. 

The flexibility of the reflector did affect the vibration of the boom and was related to 
pointing error. The multibody system allowed the effects of reflector flexibility to be 
investigated without reanalyzing the entire system. 

The program also applied easily to different orientations of the reflector. Modal anal- 


ysis did not have to be repeated as in the case of a single-body model. 


B. RECOMMENDATIONS 

Application of this research could include a study of the effects of a full range of dif- 
ferent configurations on the system dynamics. Also, damping effects can easily be incor- 
porated in this work to study the dynamic interaction of each subsystem. Extension of this 
research could include a study of the effects of a flexible body connected to the spacecraft’s 


main body on three-dimensional motion and control. 
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APPENDIX A 
MODAL ANALYSIS OF REFLECTOR AND BOOM 


Figures A.1 through A.10 show the undeformed as well as the first four mode shapes of the 


reflector and boom. The figures are on pages 84-93. 
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APPENDIX B 
REPRESENTATIVE EFFECT OF REFLECTOR FLEXIBILITY 


Figures B.1 through B.10 show representative effects of two different reflector flexibilities 


on the dynamic response of the system. The figures are on pages 96-105. 


47 


APPENDIX C 
REPRESENTATIVE EFFECT OF IN-PLANE REFLECTOR ROTATION 


Figures C.1 through C.12 show representative effects of three different in-plane orientations 


of the reflector on the dynamic response of the system. The figures are on pages 106-117. 
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APPENDIX D 
REPRESENTATIVE EFFECT OF OUT-OF-PLANE REFLECTOR ROTATION 


Figures D.1 through D.8 show representative effects of an out-of-plane orientation of the 


reflector on the dynamic response of the system. The figures are on pages 118-125. 
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APPENDIX Ee 


PROGRAMS 


id nastran Natalie 

sol 3 

time 500 

cend 

title=Nomal Mode Analysis of LFMR Boom 
subtitle=Research on N-ROSS DYNAMICS 
method=10 

displacement=all 

spc=100 

CuEDUt (plo) 

set 1 = all 

GYigin 5) mor 

axes y,Z,xX 

view -302,20-.,0r 

ptitle=Undeformed Boom 

find scale outagin o,sct 1 

plot set 1,label grid 

ptitle=Mode Shape 

maximum deformation 15 

find scale,origin 5,set l 

plot modal,deformation,0,1 THRU 10,set 1,pen 4,shape 
begin bulk 

$ 

SEIGR,sid,METHOD, f1,£2,ne,nd,,,+ 

S + NORM a. e¢ 


5 

elgr,10,/MGavy GO 100.0. 80 

$ 
SCORD2C,cid,rid,al,a2,a3,b1,b2,b3,+ 
S + Cl GZ2.¢3 

5 

GOrdZy ,1, Ose Cow Oe Oe On eet 
ter 1 Ole 

$ 
SGRID,id,ep xl x2 x3 ,,¢dy ps, seid 
: 


grid...) 1 Ome Oe on 

=,*(1),=,%*(10.1953) ,=,*(14.83427) ,== $ 
=7 

grid,10,1,91.75781,0.,133.5084 
=,*(1),=,*(-13.1085) ,=,*(12.3356),== $ 


=6 

grid, 18,10) 2.0.00. 

$ 
SCBEAM,eid,pid,ga,gb,g0/x1,x2,x3 
9 


ebeam. 1}. 3 lezen 
=,*(1),=,%(1),*(1),5= § 
-14 


90 


S 

9 

SPBEAM,pid,midl,A,1I1,12,J,nsm 

S 

peeomec i leteoo, |. 2231) laz22l), 2.4462 


SMAT1 ,mid,E,G,Nu,Rho,A,Tre£,Ge 
° 

$ A = Thermal Expan. Coef. 

$ Tref = Ref. Temp. 

$ Ge = Damping Coef. (= 2 x C/Co ) 
$ 

foe Olas, 0.25, 1522.-7 

$ 

eee e1d,¢,cid,m,*xl,*x2,%3,,+ 
cere Peel 3i- 132,133 

S 

conm2,50,17,1,9.4502-2 

S 

Seelero07123456,1,18 

$ 

enddata 


1 


id nastran Natalie 

sol 3 

time 500 

cend 

title=Nomal Mode Analysis of LFMR Reflector 
subtitle=Research on N-ROSS DYNAMICS 
method=10 

displacement=all 

spc=100 

QUEDUE(D Toe) 

set 1 = all 

Origin «5 70n7 er 

axes 2,X,Y 

View 3472) 7-2 57 nO. 
ptitle=Undeformed Reflector 

find scale,origin 5,set l 

plot set l,label grid 

ptitle=Mode Shape 

maximum deformation 15 

find scale,origin 5,set l 

plot modal,deformation,1 THRU 10,set 1,pen 4,shape 
begin bulk 

$ 

SEIGRe Sic piemnOD, flat 2 mendes: 

$ + NORIT ac 


$ 

e2qr , 1O-MGly 040, 100-0, 10 

° 

SCORDZC -€id mend: al az ta3 bl b2 bs 
S femme ?.,C3 

s) 

COTdOLn Ne O- Oe) 2 0-50. oy tz 

A Gate eee, Oe IL © 


COvd2 C62 Orr Oe sy OO 0 cee ced 
+ab, 10.70 2 23- 


$ 

SGRID ,1d cole x2 5%5 Car ps se la 
$ 

grid, l, 27 iewew >. 40es 
= *(1),=,=,*(30.),== 

=10 
Grids 72 occ Oro, Soe 
=,*(1),=,=,%*(30.),== 

=10 
Grid,25,2,60.60)5.,, 2005 
=,*(1),5,5,*(30.), = 

=10 
grid,37,2,31.5,15.,30.13 
=,*(1),=,=,*(30.),== S 
=10 
jrid,49,2755.6,15.74.72 
=,*(1),=,=,*(30.),== 

=10 


grid,61,1,0.,0.,31.67 
52 


Gamera lO. ,0.,27.65 
Gitevios,1,0.,;0..,0 
pimeyc4,1,0.50.,50. 
Gielen Oo 2,00.,0.,0. 

$ 


Seep Alyeld,p1d,ga,gb,g0/xl,x2,x3 
$ 

cbeam,1,1,1,13,64 
lees ly (1) == 5 

=10 

cbheam,13,1,13,25,64 

eee lee) == S$ 

=10 

SCAN, 2o.lrco,o) 04% 
ea. AG). FG)ia= $ 

=10 

cbeam,37,1,37,62,64 
=,*(1),=,*(1) ,== S 

=10 

cbeam,49,2,63,62,65 
ebedm, 00 ,2,02,01,095 

S 

erog,.,10,25,49 
ed by 1), == 2 

=10 
Grod 1s. 1151.49 
aes oe (1), 
=10 
emer eo 25,63 
=,*(1),=,*(1),5= $ 
=10 
Groees7 1 ).49,63 
Poe. Gh) g== S 
0 


iI 
4 
tr 


i 


emed49,11,1 2 
Oy 2 lO eee 
=9 

erod7o0 <i i2 1 
enod 01.) ieleeks., 7.4 
aol) = ele (1), 
=9 

Glod),./ 2,1 174,13 
erode 737) 125.26 
al perme (1), 
=9 

Groa,c4 1 1e35,2> 
God .co ine 36 
=,*(1),=,*(1),*(1) , == S 
=9 

erod,96,11 748,37 
GWwOadj 27 - li 25,50 
yar (1), == $ 
=9 

bemed 0S ,11,36,49 


ii 
I) 
ir 


i 
il 
4 


i 
i 
Ur 
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eGrod, 1097 11497 Ze 
= (i=, Ol) ee — S 


=9 

Erod, 120.1 oum2s 

$ 
SPBEAM,pid,midl1,A,I1,12,J,nsm 
$ 


pbean, 1,1 @21452 70 01647 (0.01647 203235 
poeam, 2,17 3-927 7571907 , 3.19077 Ges ore 
prod, LO; Oeteeeo lh .6o2-3 0 
prog,11, loge o4 ie ois 10 

$ 

SMATi,mid,E,G,Nu,Rho,A,Tref,Ge 

$ 

S A = Therniaieixpan. Coer. 

S Tref = Ret. Temp. 

S$ Ge = Damping Coef. (= 2 x C/Co ) 

$ 

Matl, 1, Wolo 22s oz e ees 

9 

SCONMZ Festa ne cree nl kc ee 

Ss +, ll) 2 ol bec ss 


S 

ConmZ 1) ez) 4502-5 
=,*(1),*(1),== 

=10 

eonm2 13713 2,070%52=5 
=,*(1),*(1),=5 

=10 

econn2, 25,2572 ,0nel22=s 
=,*(1),*(1),5= 

=10 
conmz2,37,37,2,0.4438-3 
=) =a 

=10 


conm2,49,62,1,0.626-2 
conm2,50,63,1,0.649-2 - 

5 

spel ,100, 123456763 etinuses 
5 

enddata 
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KERR KK KKK KAKA RK KK AKA AR RK KAR KKK K KARR KKK KKK KKK KARA KKK KKK KAKA KKK KK KKKKKAKAKER 


NPT,NPTC 
Me 

J3 

B2MA 

|e aang 
Me ee 8 


AHASS 
CMASS 


BMASS 
5 (2. 
DX 


Fie ates FIBaZ 
Peer ter ,F1CZ 
Poe rola, rsi3 


ONGAB , ONGAC 

Poe 2c 

PSESIN 

Bilt 

PEAY 

PHZS 

CMS 

ALSHS 

B2ous 

X3SMS 

mec 

eel 

ZC 

mieZiS 

m2 oli 

RUSS 

eee isos ane gare 
mono, © bi0,P B25 

PRAY 


tee Oe Oe OR OR OF ORR Ae FF OR eR or a oe OO A A A Oe ee ee 


BOOH AND REFLECTOR. 


SF INS AAES 


dee eee ueetO THE DISPLACEMENT OF 
BODY a vowes ER OF MASS FROM THE SPIN AXIS 
jd COORDUNATES OF THE GRID POINTS ON 


Beby B 


ene weCOORD Eales OF THE GRID POINTS ON 


BODY 1G 


THIS PROGRAM IS DESIGN TO READ IN MODAL INFORMATION FROM 
A NASITRAN DATA FILE AND CALCULATE TIME CONSTANTS FOR A DSL 
Pr See lO FILES ARE READ CONTAINING THE DATA FOR THE LFHMR 
A PARAMETER STATEMENT IS USED TO DEFINE 
THE NUMBER OF GRID POINTS AND MODES FOR BOTH BODIES IN THE 
BEGINNING FOR A MORE GENERAL PROGRAM. 
THE FOLLOWING PARAMETERS ARE USED IN THIS PROGRAM 
olde ieeemOr GRID POINTS FOR BODY B AND C 
RES oe ely 
Site etUMeeeeor ITODES FOR BODY B AND C 
REoheet fiero BE USED ITO CALCULATE 
LHe OMen tT OR SINERTIA OF BODY A ABOUT THE 


: THE TOTAL MASS OF BODY A 


: Ali ARRAY OF THE MASS AT EACH GRID POINT OF 


BODY © 


: A DUMMY ARRAY TO READ IN MODAL DATA 
MlttewocomroR ULT LENGIH OF BODY B 
es eewiiea ello tANCE BETWEEN GRID POINTS ON 


BODY B 


mo nebo Ore SCDY 8S TRANSLATIONAL MODE 
PheapeowROWS ARE DOFS, COLUMNS ARE MODES) 
Wel Chs OF sBGOye8 S TRANSLATIONAL MODE 

Si oeeom ccheeroniead AS FIBX,FIBY,FIBZ ) 

: ARRAYS OF THE ROTATIONAL MODES OF THE HINGE 
BOD ee AND BODY C 


POINT BEIVWEEN 
; [HES AUGULAR 
INTEGRAL OF 
EATEGRAL OF 
INTEGRAL OF 
INTEGRAL OF 
INTEGRAL OF 


: SUMMATION 

SUMMATION 
: SUMMATION 
: SUMMATION 
: SUMMATION 
: SUMMATION 
: SUMMATION 
: SUMMATION 
: SUMMATION 
: SUMMATION 
: SUNHNATION 


OF 
CE 
Cr 
GF 
OF 
Ole 
OF 
OF 
OF 
OF 
OF 


FREQUENCY OF BODIES B AND C 
P1**2+P2**2*RHO*DX 

(PL*FIBX + P2*FIBY)*RHO*DX 
(PIXFIBY + P2*FIBX)*RHO*DX 
(FIBX(I,J)*FIBY (I,K) *RHO*DX 
FIBZ(1,J)*FIBZ(I,K)*RHO*DX 
GRID POINT MASSES OF BODY C 
X1**2 * CMASS 

X2**2 * CMASS 

X3**2 * CMASS 


KIS SeCuASS 
NZeoeeASS 
oe CHASS 
MIRZA CMASS 
Meo SHASS 
iS cCiaASs 


Puen 2) CMASS 
ec imOlmermr Derr wga(. J) )~FICK,Y,2(1,K)*CMASS 
eOUMMATUONSOF FICM@l,J)*FICY(I,K)*CMASS 


515, 


+ 3 OF FF + FF OF Fe OF OF OM OU Um OM OU ON Oe Oe Oe Ue Ue 


+ OF OF OF OF OF OF OOOO HH 


x KPPBX yee : SUMMATION OF X1*FICX,Y,Z*CMASS x 
* Re E Key 7 : SUNMATION OF X2*FICX,Y,Z*CMASS * 
x KSPBR YZ : SUMMATION OF X3*FICX,Y,Z*CMASS * 
* Caen : A(1,1) CONSTANT TERM * 
* U1D12 : A(1,1) COEFFICIENTS OF BODY B GENERALIZED * 
* COORDINATES x 
* U1D13 : A(1,1) COEFFICIENTS OF BODY C GENERALIZED x 
* COORDINATES a 
* UJD11 : A(1,J+1) * 
* U1D21 : A(1+J,1) COEFFICIENTS OF BODY B GENERALIZED * 
x COORDINATES * 
* U1D22 : A(1+J,1) COEFFICIENTS OF BODY C GENERALIZED * 
* COORDINATES * 
* UKD11 : A(1,1+M+K) x 
* Ubieysiat : A(1+M+K,1)COEFFICIENTS OF BODY B GENERALIZED * 
* COORDINATES x 
* U1D32 : A(1+M+K,1)COEFFICIENTS OF BODY C GENERALIZED * 
x UJp2l : A(1+I,1+J) x 
* UKD21 : A(1+J,14+M+K) * 
* B12UUJ : B(1) COEFFICIENTS OF 2U1 TIMES BODY B'S * 
* GENERALIZED SPEEDS * 
* B12UUK : B(1) COEFFICIENTS OF 2*U1 TIMES BODY C'S * 
* GENERALIZED SPEEDS * 
* BU221 : B(1+J) U1**2 CONSTANT TERM x 
* BU222 : B(1+J) COEFFICIENTS OF U1**2 TIMES BODY B'S * 
* GENERALIZED COORDINATES * 
* BU223 : B(1+J) COEFFICIENTS OF U1**2 TIMES BODY c'S * 
* GENERALIZED COORDINATES * 
* : ALSO B(1+M+K) COEFFICIENTS OF U1**2 TIMES * 
* BODY B'S GENERALIZED COORDINATES * 
* B22UUJ : B(1+J) COEFFICIENTS OF 2*U1 TIMES BODY B'S * 
x GENERALIZED SPEEDS * 
* B22UUK : B(1+J) COEFFICIENTS OF 2*U1 TIMES BODY c'S * 
x GENERALIZED SPEEDS * 
* B3U211 : B(1+M+K) CONSTANT COEFFICIENTS OF U1**2 * 
* B3U213 : B(1+M+K)COEFFICIENTS OF U1**2 TIMES BODY c's * 
* GENERALIZED COORDINATES * 
* B32UUJ : B(1+M+K)COEFFICIENTS OF 2*U1 TIMES BODY B'S * 
* GENERALIZED SPEEDS * 
* B32UUK : B(1+M+K)COEFFICIENTS OF 2*U1l TIMES BODY c'S * 
* * 


GENERALIZED SPEEDS 
FPR KKK IK RK KEK IKK KKK KKK RK KERR KKK KKK KKKKKKKKEKKAKKKKKAKKKKKKKKRARARKEKE 
IMPLICIT REAL*8(A-H,0-Z) 
INTEGER NPT,NPTC,M,MC,I,J,K 
PARAMETER (NPT=17,NPTC=63 ,M=10,MC=10) 
REAL*8 J3 
DIMENSION P1(NPT),P2(NPT) ,P3(NPT),P1SQ(NPT) ,P2SQ(NPT) , RPFMIX(NPT ,M 
+), P1SP2S(NPT) ,RPSSQ(NPT) ,CMASS(NPTC) ,X1(NPTC), 
+X%2(NPTC) ,X3(NPTC) ,FIBZ(NPT,M),PHI2C(NPTC), 
+FIBX(NPT,M),FIBY(NPT,M),FICX(NPTC,MC) ,FICY(NPTC,MC) , FICZ(NPTC 
+, NC) ,PLFBX(NPT,M) ,P2N(NPT) ,P2FXN(NPT,M), PLFY(NPT,M) , PFIMIX(NPT,M) 
+, P2FBY(NPT,M),P1P2FA(NPT,M),RPSFA(NPT,M),PSFSIN(M),PSI1(M),PSI2(M) 
+,PSI3(M),FBX(MC) , PBY(MC) , X1PBX(MC) ,X2PBY(MC) ,X2PBX(MC) , X1PBY (NC), 
; +X1PBZ(MC) ,X2PBZ(MC) ,X3PBX(MC) ,X3PBY(MC) ,X3PBZ(MC), 
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+U1D12(M),ULD13(MC) ,UJD11(M) ,UKD11(MC),B12UUI(M) ,BL2UURK (NC), 
+BU221(M),B3U211(MC) ,U1D21(M,M),U1D22(M,MC),UJD21(M,M),UKD21(M,MC), 
+BU222(M,M) ,BU223(M,MC),B22UUJ(M,M),B22UUK(M,MC) ,B32UUJ(M,MC), 
+U1D31(N,MC),B3U213(NC,MC),B32UUK(MC,MC), 
+PHIPHI (NPT) ,PHEXY(M,M),PHZS(M,M),PBZ(MC),PBXS(MC,MC), 
+PBYS(MC,MC) ,PBXY(MC,MC),BIN1I(M),U1D32(MC,MC) ,PBZS(MC,MC), 
+FREQB(M) ,OMGAB(M), FREQC(MC) , OMGAC(MC) , BMASS (NPT) 
DATA J3,B,AMASS,DX,RHO/12.0,12.,.1295,18. ,1.80631D-4/ 
* READ IN CONSTANTS 
CALL MODATB(M,NPT,P1,P2,P3,FIBX,FIBY,FIBZ,PSI1,PS1I2,PS13,BMASS, 
+FREQB, ONGAB) 
CALL MODATC(MC,NPIC,X1,X2,X3,FICX,FICY,FICZ,CMASS, 
+FREOC, OMGAC) 
* CALL ALL CONSTANT SUBROUTINES 


CALL BINT11(NFT,RHO,DX,P1,P2,P1S0,P2SQ,P1SP2S,RPSSO,P1P2IN) 
CALL BINT12 (NPT,M,P1,P2,FIBX,FIBY,PIFBX,P2FBY,P1P2FA,RPSFA, 
+RHO,DX, PSFSIN) 

CALL BINT1I (NPT,M,P1,P2,FIBX,FIBY,RHO,DX,P2N,P2FXN,PI1FY, 
+PFIMIX, RBPFMIX,BINII) 

CALPXY(NPT,M,FIBX, FIBY,RHO,DX,PHIPHI, PHBXY) 
FIBZSQ(NPT,M,FIBZ,PHIPHI ,DX,RHO,PHZS) 
SUNMASS(NPTC,CMASS,CMS) 


CALL 
CALL 
CALL 
CALL 
Cae 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
Cane 
CALL 
CALL 
CALL 
Cues 
CALL 
CALL 
CALL 
CALE 
CALL 
CALL 
CAL 

CALL 
CALL 
CALL 
CALL 
BZA 


XSOCHS 
“SOCHS 
XSQCMS 


(NPTC,CMASS,X1,X1,X1SMS) 
(NPTC,CNASS ,X2,X2,X2SHS) 
(NPTC,CMASS ,X3,X3,X3SMS) 


SUNG rier ae CASS, x1Ci> 
Sues (NF Ter CASS, x2CM) 
SUNG. (NPIC RAs eHASs ,43CN ) 


ZSOCHS 
XSOCHS 
XSOCMS 
CALPBM 
CALPBi! 
CALPBIf 
PBSYMS 
PBSYMS 
PBSYMS 
CALPBS 
CALXPB 
CALXPB 
CALXPB 
CALX?B 
CALXPB 
CALXPB 
CALXPB 
CALXPB 
CALKPB 


(HET@G GMASS M1, a2, 812HS) 

(HPTC /GHASS (XZ, 43 ,XZ23H5) 

CIPI GVGMASS, x1 ,43,X13iS) 

CWRU eC BIC, CHASS PBA) 

(NEG UC EP RGYAGhiASSsrre 

ee GC ier CZICMASS , PBZ) 
ipewie Prem ielios | HiZze, PBNS) 
Cie Gye nley Ciass FHL 26, FBYS) 
Cite ieypleZ CiASs -PHEZC, PBZS) 
CNETG Ney phew nile (GMASs,PHIZC,PBXY) 
OIG Geerrne x, CMASS , XI PBX) 
CP TC icy dae iCy (CHASS -X1PBY) 
(EITC M@rah aL C2 CHASS XIPBZ) 
CNETG AH pee tex, CNASS, 22PBX) 
Ce eMC eZ ey, CHAS Se AZPBY ) 
Ci Te sie .2 fe7 CASS ,aZPBZ) 
MeO Me 53, fF ICX,/CHASS , XSPBR) 
Miele te, Koy ewer, CMASS ;-X3PB%) 
Cee ie, ney FC Z, Chass | XSPBZ) 


= B*B”AMASS 
HIH2CM = (PL(NPT)*PL(NPT) + P2(NPT)*P2(NPT))* CHS 
U1D11 =J3+B2MA+P1P2IN+H1H2CM+X1SMS+X2SNS+2.*(P1 (NPT) *X1CM+P2 (NPT) 


+ *¥2CM) 


DerlOsd =), 4 


a7 
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r2Z 


FXTER = FIBX(NPT,J)*(P1(NPT)*CMS + X1CM) 

FYTER = FIBY(NPT,J)*(P2(NPT)*CMS + K2CM) 

SIITER = -PSI1(J)*(P2(NFT)*X3CM + X23MS) 

SI2TER = PSI2(J)*(P1(NPT)*X3CM + X13NS) 

SI3TER = PSI3(J)*(P2(NPT)*XICM - P1(NPT)*X2CM) 
BU221(J) = PSFSIN(J)+FXTER+FYTERtSI2ZTER+SI3TER+SIITER 
UIDI2Z(J) = 2. = eucennen 

CONTINUE 


DO 11K =1,mMc 
U1D13(K)=2.*(P1(NPT)*PBX(K) +X1PBX(K) + P2(NPT)*PBY(K) + X2PBY(K)) 
CONTINUE 


DO 12 J =1,M 
FXTER = FIBX(NPT,J)*(P2(NPT)*CMS + X2CM) 


FYTER = FIBY(NPT,J)*(PL(NPT)*CMS + X1CM) 

SIITER = PSI1(J)*(P1(NPT)*X3CM + X13MS) 

SI2TER = PSI2(J)*(P2(NPT)*X3CM + X23MS) 

SI3TER = PSI3(J)*(P1 (NPT) *X1CM+K1SMS+X2SMS+P2 (NPT) *X2CM) 


UJD11(J) = BIN1I(J)-FXTER+FYTER-SI1TER-SI2TER+SI3TER 
CONEINUE 


XCALCULATES A(1,3) AND PART OF A(3,1) 


i 


14 


LS 
A 1l+J 


ry 
16 


pO 13 k= 1c 

UKD11(K) = -PBX(K)*P2(NPT)-X2PBX(K)+PBY(K)*P1(NPT) + X1PBY(K) 
CONTINUE 
ere 4 SS 

FXTER = FIBX(NPT,J)*(P1(NPT)*CMS + X1CM) 


FYTER = FIBY(NPT,J)*(P2Z(NPT)*CMS + X2CM) 
SIITER = PSI1(J)*(P2(NPT)*X3CM + X23MS) 
SIZTER = PSI2(J)*(P1(NPT)*X3CM + X13MS) 
SISTER = PSI3(J)*(PL(NPT)*X2CM + P2(NPT)*K1CH) 


B12UUJ(J)=-PSFSIN(J)-FXTER-FYTER-SI1TER+SI2TER-SI3TER 
CONTINUE 
DONS sh eG 
B12UUK(K)=-P1 (NPT) *PBX(K)-X1PBX(K)-P2 (NPT) *PBY(K)-X2PBY(K) 
B3U211(K) = -B12UUK(K) 
CONTINUE 
EQUATION 


DO) heat ale 
DOI J=1 ii 


FXTER=FIBX(NPT,J)*(PSI1(1)*X3CM-PSI3(1)*X1CM-FIBY(NPT,I)*CHS) 
FYTER=FIBY (NPT ,J)*(PSI2(1) *X3CM-PSI3 (I) *X2CM+FIBX (NPT ,I) *CHS) 
SIITER=PSI1(J)*( (PSI3(I)*P2(NPT) + FIBX(NPT,I))*X3CM + 


+(PSI1(I)*P1(NPT) + PSI2(I)*P2(NPT))*X2CM + PSI2(1)*(X3SMS+X2SMS) 
++PSI1(1I)*X12NS) 


SI2TER=PSI2(J)*( (PSI3(I)*P1(NPT) - FIBY(NPT,I))*X3CM + 


+(PSI1(I)*P1(NPT) + PSI2(I)*P2(NPT))*X1CM + PSI1(I)*(X3SMS+X1SMS) 
++PSI2(I)*X12MS) 


SISTER=PSI3(J)*( (-PSI3(I)*P1(NPT) + FIBY(NPT,1I))*X2CM + 


+(PSI3(I)*P2(NPT)+FIBX(NPT,1I))*X1CM-PSI1(1I)*X23MS+PSI2(1I)*X13MS) 


U1lD21(J,1)=PHBXY(1,J)-PHBXY (J, I )+FX1ER+PYIER@ Sl ee oe Poet oe ee 
CONTINUE 
CONTINUE 


08 


12 
18 


Zu 
20 


Doms gra it 

DO 19 K = 1,MC 
FXTER=FIBX(NPT,J)*PBY(K) 
FYTER=FIBY(NPT,J)*PBX(K) 
SIITER=-PSI1(J)*(X3PBX(K)+P1 (NPT) *PBZ(K)+X1PBZ(K) ) 
SI2TER=-PSI2(J)*(X3PBY(K)+P2 (NPT) *PBZ(K)+X2PBZ(K) ) 
SI3TER=PS13(J)*(2.*(X2PBY(K)+XLPBX(K) )+P2(NPT)*PBY(K)+P1(NPT)* 


- PBX(K)) 
U1D22(J,K)=-FSXTER + FYTER + SIITER + SI2TER + SI3TER 
CONTINUE 

CONTINUE 


DO 20 i= 1,M 
BC Di aa 
FXTER=FIBX(NPT,J)*(PSI2(I)*X3CM - PSI3(I)*X2CM) 


cf Pete Nei) Ceo) seins I 303 )*XZ2CM) 


FYTER=FIBY(NPT,J)*(-PSI1(1)*X3CM + PSI3(I)*X1CM) 


+ +FIBY(NPT,1I)*(-PSI1(J)*X3CM + PSI3(J)*X1CM) 


FZTER=FIBZ(NPT,J)*(PSI1(1)*X2CM - PSI2(I)*X1CM) 


+ +FIBZ(NPT,I)*(PSI1(J)*X2CM - PSI2(J)*X1CM) 


Sele Rese (sli 14 xeSiStaAZsMs )=PSI2(1)"x12Ms) 


+ +PSI1(1)*(PSI1(J)*(X3SMS+X2SMS ) -PSI2(J)*X12MS) 


SI2TER=PS1I2(J)*(PSI2(1)*(X3SMS+X1SMS ) -PSI3(1I)*X23MS) 


n +PSI2(1)*(PSI2(J)*(X3SMS+X1SNS ) -PSI3(J)*X23MS) 


SI3TER=PS1I3(J)*(PSI3(I)*(X2SNS+X1SMS )-PSI1(1I)*X13HS) 


te +PSI3(1I)*(PSI3(J)*(K2SNS+K1SHMS )-PSI1(J)*xX13MS) 


IF(I .EQ. J)THEN 
UJD21(1I,J)=1.+FXTER+FYTER+FZTER+SIITER+SI2TER+SI3TER 

ELSE 
UJD21(1,J)=FXTER+FYTER+FZTER+SI1TER+SI2TER+SI3TER 
UJD21(J,1)=UJD21(I,J) 

ENDIF 

CONTINUE 

CONTINUE 


SCALCUMEMES A(Z,3) AND A(3,2Z) 


Wee Kk = 1 Ne 

DO 23 J = 1,M 
FXTER=FIBX(NPT,J)*PBX(K) 
FYTER=FIBY (NPT ,J)*PBY(K) 
FZTER=FIBZ(NPT,J)*PBZ(K) 
SIITER=PSI1(J)*(X2PBZ(K)-X3PBY(K) ) 
SIZTER=PS12(J)*(X3PBX(K)-X1PBZ(K) ) 
SI3TER=PSI3(J)*(X1PBY (K)-X2PBX(K)) 
UKD21(J,K)= FKTER+FYTER+FZTER+SIITER+SI2TER+SI3TER 
CONTINUE 

CONTINUE 


Domes J 20a 
en25.1 = Ion 
FATER=FIBX(NPT , J)*(-PSI3(1)*X2CM+PS12(1)*xX3CM) 
FYTER=FIBY (NPT ,J)*(-PSI1(1)*X3CM+PSI3(1I)*X1CM) 
SIITER=PSI1(J)*(-FIBY(NPT ,1)*#%3CM+PSI1(1I)*(X3SMS-X2SMS-P2 (NPT) 
+*X2CN)+PS12(1)*(P1 (NPT) *X2CM+K12NS )+PSI3(1)*P1 (NPT )*X3CM) 
SIZTER=PS1I2(J)*(FIBX(NPT,1I)*X3CM+PS1I2(1)*(X3SMS-X1SMS-P1 (NPT) 


og 


+*X1CH)+PSI1(1)*(P2(NPT)*XLCM+X12NS )+PSI3(1)*P2(NPT)*X3CM) 
SI3TER=PS13(J)*(-FIBA(NPT,1I) *X2CN-PS13(1)*(P2(NPT)*X2CM+P1 (NPT) 
+*K1CM)-PSI1(1)*(X13MS ) -PSI2(1)*X23MS+FIBY(NPT,1I)*X1CM) 
IF (I .EQ. J) THEN 
BU222(J,1)=1.-FPHZS(J,1)-CMS*FIBZ(NPT,1)*FIBZ(NPT,J)+FXTER+FYTER 


+ +SILTER+SI2TER+SI3TER 
ELSE 
BU222(J,1)=-PHZS(J,1)-CMS*FIBZ(NPT,1)*FIBZ(NPT,J)+FXTER+FYTER 
s: +SI1TER+SI2TER+SI3TER 
ENDIF 
25 CONTINUE 


24 CONTINUE 


DO 26 K=1,MC 
DO 27s 
FXTER=FIBX(NPT,J)*PBX(K) 
FYTER=FIBY(NPT,J)*PBY(K) 
SI1TER=-PSI1(J)*(X3PBY(K)+P2(NPT)*PBZ(K)+X2PBZ(K) ) 
SI2TER=PSI2(J)*(X3PBX(K)+P1 (NPT) *PBZ(K)+X1PBZ(K) ) 
SI3TER=PSI3(J)*(P2(NPT)*PBX(K)-P1(NPT)*PBY(K) ) 
BU223(J,K)=FXTER+FYTER+SI 1TER+SI 2TER+SI3TER 
1) CONTINUE 
26 CONTINUE 


DO! 2ema et 

DO 29 I=1,M 
FXTER=FIBX(NPT,J)*(FIBY(NPT,1I)*CMS+PS13(1)*X1CM-PSI1(1)*X3CM) 
FYTER=-FIBY(NPT,J)*(FIBX (NPT ,1I)*CNS+PS12(I)*X3CM-PSI3(1I)*X2CM) 
FZTER=FIBZ(NPT,J)*(.5*PSI3(1)*X3CM) 

SIITER=PSI1(J)*(FIBX(NPT,1)*X3CM+PSI2(1)*X3SMS-PSI3(1)*1.5*X23MS) 

SIZTER=PS12(J)*(FIBY(NPT,1)*X3CH-PSI1(1)*X3SMS+PS1I3(I)*1.5*X13MS) 

SI3TER=PS13(J)*(-FIBY(NPT,1I)*X2CM-PSI2(1)*X13MS+PSI1(1)*X23MS- 


i: FIBX(NPT,1I)*X1CM) 
B22UUJ(J,1)=PHBXY(1I,J)-FHBXY(J,I) + FXTER + FYTER + FZTER + 
ie SILTER+SI2TER+SI3TER 
29 CONTINUE 


28 CONTINUE 


DO 30 K = 1,MC 
pO 31 J=1,M 
FATER= FIBX(NPT,J)*PBY(K) 
FYTER= -FIBY(NPT,J)*PBX(K) 
SIITER= PSI1(J)*X3PBX(K) 
SIITRK= SILTER -PSI1(J)*(P1(NPT)*PBZ(K)+X1PBZ(K) ) 
SI2TER= PSI2(J)*X3PBY(K) 
SI2TRK= SI2TER -PSI2(J)*(P2(NPT) *PBZ(K)+X2PBZ(K) ) 
SI3TER= -PSI3(J)*(X1PBX(K)+X2PBY(K) ) 
SI3TRK = PSI3(J)*(P1(NPT)*PBX(K)+P2(NPT)*PBY(K)) 
B22UUK(J,K) = FXTER + FYTER + SIITER + SI2TER + SI3TER 
B32UUJ(J,K) =-B22UUK(J,K)+.5*%3PBZ(K)*PSI3(J) 
U1D31(J,K) = FXTER + FYTER + SIITRK + SI2TRK + SI3TRK 
31 CONTINUE 
30 CONTINUE 
*  14M+K EQUATION (UKD11 IS 1RST TERM; U1D31 IS 2ND;IDENTITY IS 3RD) 
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x B63 (B3U211 IS IRST TERM;BU223 IS U2 Q(I) TERH; 


* B32UUJ IS 2U *U(1+1)TERM 
Doms2zek = 1Ne 
DO 33 J =K,NC 

IF (K .EQ.J) THEN 


L 


= PBZS(J_,K) 


B3U213 (J,K) = 
ELSE 
B3U213 (J,K) = -PBZS(J,k) 
Bauz13(K,d) = Bsuzis(s,K) 
ENDIF 
33 CONTINUE 


a2 CONTINUE 
DO sae Ne 
DOgso DO = tC 


B32UUK(J,K) = PBXY(K,J)-PBXY(J,K) 
UIDS2(3,K) = - B32UUK(J,K) 


55 CONTINUE 
34 CONTINUE 
Ceo Vee VieeRE SULTS IN FILE Z5 


em (25, '(10K,*'ALIC'*)") 


WRITE (25,100) U1D11 
100 FORMAT (D15.6) 


Wemmece> (107%, “VALIO"" 9X, “MATIC'’ 8x,''B1QO"™ 9x, ''BUJUILC'’ , 9X, 


+''OHGAB'')') 
DO 1 J= 1," 


WRITE(25,200)J,U1D12(J) ,UJD11(J) ,Bl2UUd (J) ,BU221(J) , OMGAB(J) 


1 CONTINUE 
200 FORNAT (14,5D13.6) 


Pee eee elon BOB tyex S ALIOBTY Ox, 9 'ALKC''9X,''BRKULC'' ,9X, 


tee OGAC* * )) ) 
BOeZz k= 1,0C 


Eee 25, 200) K BIZUUK(K) (UlD13( Kh) UKDI1(K),B3U211(K) ,OMGAC({K) 


Z CONTINUE 


mene ole RIC’! lik eo een,  BIULO'' ,13X,'"BJ2UQ'')') 


BOwo. t= 1K 
D0 4JJ=1,i 


Heme s)500)),1,UJ021() 1) UlD21(J 1) ,BuUz22(J,1I) ,B22UUJ(J,1) 


4 CONT EIVUE 
3 CONTINUE 
500 PORIIEa 224 ,4D15.6) 


Piezo (tzn,' AKIOB'!', 9X%,''BKULOB'' ,9X,''BK2UQB'')"') 


DOM Ka IC 
DOSS e= 1 iC 


hoe 200)1,K,ULD32(1,K),B3UZ213(1,K),B32UUK(1,K) 


6 CONTINUE 
5 COMMIHUE 
500 PORMAT( 214..3D15.6) 


Pee leeereer eh On eon) TAKIIC!! 12K.) BKUIO'',11X,''BUZUQB' ' ) 


+!) 
DO 7 K =1,HC 
ee m eee 


Prem 257 500)d7 8, UlD22() /K) UKD2Z1(0 ,K) ,BU223(J,K) ,B22UUK(J,K) 


8 CONTINUE 
7 CONPINveE 
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Shs 
98 
400 


Sif 
96 


8 


a2 
900 


WRITE(25,'(12X,''AK1Q'' ,9X,''BK2UQ'')") 
DOWSe Ki sieie 
pO 99 J = 1,m 
WRITE (25,400) J,K,U1D31(J,K) ,B32UUJ(J,K) 
CONTINUE 
CONTINUE 
FORMAT (214,2D15.6) 
DO 94 J = 1\M 
DO 95 I =1,NPT 
WRITE(25,900)FIBX(I,J),FIBY(1I,J),FIBZ(I,J) 
CONTINUE 
CONTINUE 
DO 96 J = 1,MC 
DO 97 I = 1,NPTC 
WRITE(25,900) FICX(I,J),FICY(I,J),FICZ(I,J) 
CONTINUE 
CONTINUE 
DO 93 J = 1,M 
WRITE (25,900) PSI1(J),PS1I2(J),PSI3(J) 
CCNTINUE 
pe 2 ha aku 
WRITE (25,900) X1(I),X2(1I),X3(1) 
CONTINUE 
FORMAT (3D20.6) 
END 


KRAAAKRAKAKKKKKKAKRAKRAKARAKRKAKKEKRKKKKAKARKAKKEKKKRKKARKKAKKKKKKAKKAKRAAKAKKKKRKKAKK 


= 
x 
* 
* 


* 
VECADD x 


( ADD TWO ARRAYS TO PRODUCE AN ARRAY ) - 
k 


RA AAARARAA RAK AKAAAKARAA KARA AA AAA ARK RAK KARR KAK RRA KAKA AARKAAARKAARAKAARKAKAKK KAR 


10 


* 


SUBROUTINE VECADD (NPT,A,B,ARG) 
REAL*8 A(NPT),B(NPT) ,ARG(NPT) 
be) 0) i STIR ae 

ARG(I) = A(I) + B(I) 


CONTINUE 
RETURN 
END 
KAAEARAEKARAKAKKKEKRKKKKAKAKAKKAARKAAKKARAKAAAARKAAKKKRKEKKAARKKKAKKRKRKEKKAKKRKKARKK 
MULTSV x 
( MULTIPLIES A SCALAR TIMES A VECTOR ) - 


* 


AAARREKKKEKAKERAKKAKKAKKRARRAKKEARKKKKKARKKARKAKRRARKKKRKARKARKEKAKRKKAKAKKAKKRRRAAKKKR 


10 


* 
* 


SUBROUTINE MULTSV (NPT,A,B,ARG) 
REAL*8 A,B(NPT),ARG(NPT) 
DO 10 I = 1,NPT 

ARGO] = Ane eae 


CONTINUE 

RETURPII 

END 
RAKKRAAKARAAAKAKAAKRAKA AA KAKA RAK RR RA KARA AKKAAKAAAKKAAKARRAKRRAARKARKAKRAKAKAKR 

MATADD ; 

( ADDS TWO MATRICES TO PRODUCE A MATRIX ) : 
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Mee here Ae A AR APE ARARKAKRAKARAAKAAARKAAKRKKKRKKKKKRKARKREKRKRKKARKKKK EK 


SUBROUTINE HATADD(NPT,H,A,B,ARG) 
Pewee Ss ACIIPY 1) S¢NPT 11), ARGCNPT ,M) 
DO 10 J=1,/M 
DOezZ0 i = 1,NPT 
Poe td )os= ieee (Ld) 


20 COUT TNE 
10 CONTINUE 

RETURN 

END 
KAA RAR RRARK RK AR ARA AAR RA RKK AKA EAA RAR AR AKAAKARAARAKKKAAKRAK AAR RK AAKRKKKR KK 
3 SUBROUTINE MULTZS "3 
. ( Gabeuteteo A COLN TIHES A ROW VECTOR ) zs 


SUBROUTINE HULTZS { NPT,A,B,ARG) 
REAL*8 A(NPT),B(NPT) ,ARG(NPT) 
PO ee = NPT 
HEeG( i= ACT) ~ BCI) 
10 Coli erie 


RETURN 

END 
RARKARKAAARAAAARARKARRARKRERKRRKKKAKKEKKKARRARARARKRRKKKKKARKRRRKKKKRRERKKKEKKERK 
% MULTSI . 
x ( MULTIPLIES A SCALAR TIMES A MATRIX ) x 


AARKRRAKAKRARRRARKAAKRKARARAKARKARKAAKRKKKAKKKRERKAKKKKKKKAKKRKRKKKAKKRRERRRKRKRRKKE 


SUBROUTINE MULTSM (NPT,M,A,B,ARG) 
REAL*8 A,B(NPT,M),ARG(NPT,M) 
Porro d = 1.1 
DO 20 I = 1,NPT 
ARG(I,J) = A* B(I,J) 


20 CONTINUE 
10 CONTINUE 

RETURN 

END 
REAR AKREKRAK RK REE KKK KKK KA KK RRR KKK RK KKK KKK KR KKK KKKKKKKAKRKEKKKKKRKKRKRKKKK 
* SUBROUTINE MULTVM ‘ 
* ( MULTIPLIES VECTOR AND MATRIX TO PRODUCE A MATRIX ) * 


RAAT NAR KR RA AKARAKKAAARKKARAARKEKARKRAKRRAKRKARKAKRKKKRAKAAKKRKKAKKKKE 


SUBROUTINE MULTVM (NPT,M,A,B,ARG) 
REAL*8 A(NPT),B(NPT,M) ,ARG(NPT,I) 
Domine = 1 NET 

DO 20 J =1,N 

PeGo@)) = ACI) * B(I,J) 


20 CONTINUE 
10 CONTINUE 

RETURN 

END 
Re ee AAA A AAA ARR AR ARARKAKKRAKKKAKRARKARKAKKKRKRRKKEAAERKK 
x SUBROUTINE FIBZSQ . 
* ( CALCULATES THE INTEGRAL OF PIBZ(I,J)*PIBZ(1,K)*RHO*DX ) . 


REARAN MM RA PR ARR KKK RRR AR AP RK AA RRAR ARK RAK KKK RAKANARKKAAKKKAKKREKKKERKAKKE 


SUER oO ie SZASOGIEI MN, kieZ PHIPHI ,DxA,RHOMPHZS) 
Paes FIBA (MPD ieee (NPP ye PHZS(M iM), RHO, DX 
DOR Oe I 
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DO 20 J = KM 
DOeS0 7 i= 17NET 
PHIPHI(I) = FIBZ(I,J)*FIBZ(I,K)*RHO 
30 CONTINUE 
CALL TRAPZ(NPT,DX,PHIPHI,SUM) 
PHZS(J,K) = SUM 
PHZS(K,J)=PHZS(J,K) 


20 CONTINUE 
10 CONTINUE 

RETURN 

END 
RAEKKRKKKKRAKKKRAKRKKRKRRREKKRKRRKKRRKKRRKKRKRRRREKRKRKRRKERRKRKAKRKRKRKRKKKRKRRRRKRRKRKRRKKKRREK 
x TRAPZ . 
x ( PERFORMS INTEGRATION APPROXIMATION USING TRAPEZOIDAL RULE) * 


KRAKKKKRKKRKKEKKKKRKKKRKRKREKKKKRKRKKKARKRKRKKKRRKRKRKKKRKRKRKERKEKRKRRKRKRRKRKRKRKRKKKEKKRRERK 


SUBROUTINE TRAPZ(NPT,DX,F, INT) 
REAL*8 DX,F(NPT),INT,SUM 
SUM = 0. 
DO 10 I= 1,8 
IF(T SEO. 1 .OR2 I EOM NPD) one 


FACTOR = l. 
Bok 

FACTOR = 2. 
1a ABB 


SUM = SUM + FACTOR*F(I) 
10 CONTINUE 
INT = DKX/2. * SUM 


RETURN 

END 
RERAKKRKRAKKRRAKRRKKRKKRKAEKRKKRRKAKKRKRKRRRRKKKRKRKRKRKRKKKRRKAKRRRKRKRKRKKRRKRKKKRKRKRKKKRKAETKER 
x MATTRP x 
x (INTEGRATION APPROXIMATION OF A MATRIX USING TRAPEZOIDAL RULE) * 


KARKAKKKKKKAKRKRKRRAKKKKKKRAKKRKARKREKRRKERKKKRRKKKKKKRRKRRKRRRKRKRKRERKRRRKKKRRRRK 


SUBROUTINE MATTRP(NPT,M,DX,F,INT) 
REAL*8 DX,F(NPT,M) ,INT(MH),SUM 
Bee) i) 2 
SUM = 0. 
DO 10 I= 1,8 
IF(I .£O. 1 .OR. I .FOmes uanuel 


FACTOR = 1. 
BLSE 

PRCTOR = 72. 
BHDTE 


SUM = SUM + FACTOR*F(I,J) 
10 CONTINUE 
IND (J) = Deen 


5 COM ahi: 

RETURN 

END 
KAAK AK AAKKRKAAKAAKAKAAKAKAAAAKAK AAA AAA AAAAKEAKRARRRRERRRERERAERKAKAKKKKKKK 
a SUHC am 
ce ( PERFORMS SUMMATION OF GRID POINT MASSES Cr geC lsc UeliEn ‘ 
a AN ARBITRARY VECTOR OF EQUAL DIMENSION ) . 


ARKEKKRARKKEKKREKRKRARRRRRAKR ARERR KARR KRAKRAKR AR ARARRARKRRAN AR RATAN A  ee 


64 


SUBROUTINE SUNC(NPTC,F,CMASS, SUI) 
REAL*8 F(NPTC),CMASS(NPTC) ,SUII 
SUM = 0. 
Demonia = 1 Nic 
SUM = SUI! + F(I)*CMASS(I) 
10 CONTINUE 


RETURN 

END 
KEAK RAK AA AKA AKAA AKA AKAKAKAKK AAA A AK AKAAKAKK KARA KA KKAKKAKAKAKKAKKK KKK A KK AK RAK 
z MATSUI! = 
a ( PERFORMS SUMMATION OF THE GRID POINT MASES OF BODY C TIMES = 
a AN ARBITRARY MATRIX OF EQUAL DOF) a 


RARKKKRRKEKAKARRKEKKKRKKAAKKAKAKRAKKARRKRKKRKKKKKAKKKKKKKKKRKRKKREKRKKKKKRAKERAKK 


SUBROUTINE MATSUM (NPT,M,F,CMASS,SUM) 
REAL*8 F(NPT,M),CMASS(NFT) ,SUM(M) 
DO 10 J= 1,H 
Der 20 1 = 1, NPT 
SUM(J) = SUM(J) + F(I,J)*CMASS(I) 


20 CONTINUE 
10 CONTINUE 

RETURN 

END 
KARA KAR RK RA A KRK KARR RAK RK KKK KKRA KR KR K KK RK KARR KKK RKKRKKKRKKKK KKK KKK KKKKKAKKKK 
* BINT11 x 
* ( PERFORMS INTEGRATION OF (P1**2 + P2**2*RHO*DX) x 


KRAKKRKKKRRRRKRARKRKRRKKEKKKKRKKRKKRKRRRKKKKKKRKKKKRRAKKAKRRRKKKKKKRKRKKKKKKRKAKKKK 


SUBROUTINE BINT11(NPT,RHO,DX,P1,P2,P1SQ ,P2SQ,P1SP2S,RPSSQ,P1P2IN) 
REAL*8 P1(NPT),P2(NPT) ,P1P2IN,P1SQ(NPT) ,P2SQ(NPT) ,RPSSQ(NPT) 
REAL*8 P1SP2S(NPT) , BHO, DX 

CALL MULT2S(NPT,P1,P1,P1SQ) 

CALL MULT2S(NPT,P2,P2,P2SQ) 

CALL VECADD(NPT,P1SQ,P2SQ,P1SP2S) 

CALL MULTSV(NPT,RHO,P1SP2S,RPSSQ) 

CALL TRAPZ(NPT,DX,RPSSQ,P1P2IN) 


* 
* 


RETURN 

END 
KRKEEKARKKARKKAKKKRKKKKKKKKKKKKKKAKKKKKERKKKKKKKKKAKKKKKKKAKKKKKKKKKKAKKAKKAKRAKKKKE 
. BINT12 
x ( PERFORMS INTEGRATION OF (P1*FIBKX + P2*FIBY)*RHO*DX ) 


KARA KRAK A RAKARAKAKKKAR KAR AK AAAAKRARERKAKAR AK AKRRKK RK AKKKAKKAKAARAAKKKK KAA 
See OUieNe BINTI2 (NPT M, Pl P2 PIBX FIBY ,PIFBX,P2FBY ,PIP2FA,RPSFA, 


+RHO,DX, PSFSIN) 

REAL*8 P1(NPT),P2(NPT),FIBX(NPT,M),FIBY(NPT,M),P1FBX(NPT,M), 
+P2FBY(NPT,H),P1P2FA(NPT,M),PSFSIN(M),RPSFA(NPT,M),RHO,DX 
CALL MULTVM(NPT,P1,FIBX,P1FBX) 

CALL MULTVM(NPT,P2,FIBY,P2FBY) 

CALL MATADD(NPT,P1FBX,P2FBY,P1P2FA) 

CALL MULTSM(NPT,M,RHO,P1P2FA,RPSFA) 

CALL MATTRP(NPT,M,DX,RPSFA,PSFSIN) 


* 


RELUR 

END 
KRRERKAEKKKAKAKRKRAKKRARA RAK RAR KK AK RK KKK KK AKA KAA AAKAKRKARAKRKR KK RAK RK KAKA 
a Bl ain 
23 ( PERFORMS INTEGRATION OF (P1*FIBY - P2*FIBX)*RHO*DX ) 
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* 


KKK KRE KAKA KKK KKKR RRR RRR KKK RRR KR RRR KKK KK RR KERR RRR KK KKK RK RKRAKKKKAKAAKKHK 
SUBROUTINE BINT1I (NPT,M,P1,P2,FIBX,FIBY,RHO,DX,P2N,P2FXN,PILFY, 
+PFIMIX,RPFMIX,BINII) 
REAL*8 P1(NPT),P2(NPT),FIBX(NPT,M),FIBY(NPT,M),P2N(NPT) ,BIN1I(M) 
+, P2FXN(NPT,M),PIFY(NPT,M) ,PFIMIX(NPT,M) ,RPFMIX(NPT,M) ,RHO,DX 
Be ie) f= nea 
P2N(I) = -P2(I) 
10 CONTINUE 
CALL MULTVM (NPT,M,P2N,FIBX,P2FXN) 
CALL MULTVM (NPT,M,P1,FIBY,PI1FY) 
CALL MATADD (NPT,M,P2EXN,P1FY,PFIMIX) 
CALL NULTSM (NPT,M,RHO,PFIMIX,RPFMIX) 
CALL MATTRP (NPT,M,DX,RPEMIX,BINII) 


RETURN 

END 
KAKKKKKKKKKRKKKKKKKEKRKKKKKRKRKKKKKRKRKKKKKRKKKRKKKRRKRKKKKKRRKKKKKRKKRKKKKKKKKKK 
* CALPXY * 
x INTEGRAL OF ( FIBX(I,J)*FIBY(I,K)*RHO*DX) . 


KAKKKKKKKKRKKKRRKKKKKRKRKKKKRAKRKRRRRARKKRKRRKAKKKKRKRKRRRRRRRRRRKRRKKKKRRKKERKKR:K 


SUBROUTINE CALPXY ( NPT,M,FIBX,FIBY,RHO,DX,PHIFPHI , PHBXY) 
REAL*8 FIBX(NPT,M),FIBY(NPT,M),PHBXY(M,M) ,PHIPHI(NPT) ,RHO,SUM 
+ DA 
DO 10 K = 1,M 
DOPZ0s = 1M 
DO 30 1 = PNET 
PHIPHI(I) = FIBX(I,J)*FIBY(I,K)*RHO 
S10) CONTINUE 
CALL TRAPZ(NPT,DX,PHIPHI,SUM) 
PHB) = "SUM 


20 CONTINUE 
10 CONTINUE 

RETURN 

END 
KEKKKARARKKEKKKKKKKKKRKKKAKKKRKARKKKKKKKERRRKKKKKKKRKRKRKEKRKKKKRKKRKEKKKKKKKKKKKKEKK 
E SUMASS “ 
‘ (SUMMATION OF THE GRID POINT MASSES OF BODY C ) a 


KKKKKKKKKKKAKKAKKKKKARKKKKKRKKKRKKKKRRERKKKKEKKEKKKKKKAKRKKKKKKRKKRKKAKKKKKK 


SUBROUTINE SUMASS (NPTC,CMASS,CHS) 
REAL*8 CNASS(NPTC) ,CMS 
CMS = 0. 
DO 10h ieee TC 
CNS = CMS + CMASS(I) 
10 CONTINUE 


RETURN 

END 
KKAKKKEKRRKKKKKKEKRKKKRKKKKRKEKKKRKKAKRRKKKKRRKRKKRKKKKKRRRKKKKRRKKRRKKRKRKRKRKKKKRK 
* XSOCMS * 
x ( SUMMATION OF THE MULTIPLICATION OF TWO ARBITRARY X'S TIMES * 
* THE GRID POINT MASSES OF BODY C ) s 


KEKE KE AK KR KR RRR EK RRR KKK RRR KKK KKK RRA KKK RR KR KK KRKKEKKRKKRKKRKKRKRKRREKRKKR 
SUBROUTINE XSQCMS (NPTC,CMASS,XA,XB,XSMS) 
REAL*8 XA(NPTC) ,XB(NPTC) ,CMASS(NPTC) , XSMS 
ysMG = 10m 
BO Vout =i ane 
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Pollow=—enoio + XAG) “xBCT)“CMASS (1) 


10 CONTINUE 

RETURN 

END 
KAR KKKKEKRKKKRRARKKKKKRRRKKRKKRKRKKAKRKRKKRKKKRRRKRKRKRKKKRRKKKRKRRKRKKKKKRKKRKRRAKKKKK 
x CALPBI . 
x (SUMMATION OF AN ARBITRARY MODE SHAPE MATRIX TIMES THE GRID x 
x POINT MASSES OF BODY C ) 7 


RAKARRAARKARARRERKKKKKKRARARKKRKKRAKRRKERKRARKKRKRRKRRKRRRRKKKRKRKRKRKKARRKRKEKKKKKKKK 


SUECOUDMIESCALPBM (NPIC(HC,FIC,CMASS ,PB) 
Penta oe PGi iC MO) CuAsS(NeTe), PR(MC),SUM 
DO 1LOeerF= 1 LIC 
SUM = Q. 
DOV 2051 = 1L,NPIC 
Sure ouiet  FiGQiLaJ)-CMASS()T) 


20 CONTINUE 
PB(J) = SUM 

10 CONTINUE 

RETURII 

END 
KAKA RERARKKKRKKRKKKKKKKRKKKKREREKKKRKAREKKERKKKRKRKRRKKKRKEKKRKKKRKKKRERKKRKKKKE:K 
. CALPBS . 
k ( SUMMATION OF THE MULTIPLICATION OF TWO ARBITRARY MODE SHAPES * 
x TIMES THE GRID POINT MASSES OF BODY C ) - 


RRR RRR RAMANA KARA AKA RR ARK AKA KAR ARAMARK ARRAKRK AR AKA AKAKKKKKKKKKKEAKRRKKRRESR 


SUBROUTINE CALPBS (NPTC,MC,FICA,FICB,CMASS,PHI2C, PBS) 
REAL*8 FICA(NPTC,MC),FICB(NPTC,MC) ,CMASS(NPTC) ,PBS(MC,MC) 
REAL*8 PHI2C(NPTC),SUM 
DO 10 K = 1,MC 
POm2e 5. = lie 
DO 30 1 = irre 
PHI2C(I) = FICA(1,J)*FICB(I,K) 
30 CONTINUE 
CALL SUNC(NPTC,PHI2C,CMASS,SUM) 
FBS (J,K) = SUM 


20 CONTINUE 
10 CONTINUE 

RETURN 

END 
KAAARKKKEKARAAKRKRKKKRKEKRKRRKRKKERKKKKRKKKAKKKRKEKKAKKKKRKAKRKRKRRRKRRRREKKKKKKKEKKKK 
x PBSYMS . 
i (SUMMATION A MODE SHAPE MATRIX TRANSPOSED TIMES THE MODE SHAPE * 


e MATRIX TIMES THE GRID POINT MASSES OF BODY C)- A SYMMETRIC MATRIX * 


SUBROUTINE PBSYMS (NPTC,MC,FICA,CMASS,PHI2C, PBS) 
REAL*8 FICA(NPTC,MC),CHNASS(NPTC) ,PBS(MC,MC) 
REAL*8 PHI2C(NPTC),SUM 
DO 10 K = 1,MC 
DO 20 J = K,Mc 
DO 30 I = 1,NPTC 
Pree (l) = FICA(], J)*FICA(1,K) 
30 CONTINUE 
CALL SUMC(NPTC,PHI2C,CMASS, SUM) 
PBS (J,K) = SUM 
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FBS (K,J) =()BBomGe) 


20 CONTINUE 
10 CONTINUE 

RETURN 

END 
KRKAEKKKAKRKEKRAKKRKEKKKKKKRKRRRRRKRRKRKKKEKKKKKKKKKKKKKKKKKERKRKKAKRRKEKKKKKKKKKK 
x CALXFB . 
x ( SUMMATION OF ANY X*ANY MODE SHAPE MATRIX*THE GRID POINT MASSES * 
x OF BODY C ) x 


KAKKKKARRKRKKKKKEKKKAKKKKKKAKRRKKKRKRKRRKRRKRKKKKKEKKREKKKKRRKKKKKKKKKKKKKEKKKKKKKKKK 


SUBROUTINE CALXPB (NPTC,MC,X, FIC, CMASS,XPB) 
REAL*8 FIC(NPTC,MC) ,CMASS(NPTC) ,XPB(MC) , SUM, X(NPTC) 
Do 10 J = 1,Mc 
SUM = 0. 
DO 20 I = 1,NPTC 
SUM = X(I)*FIC(I,J)*CMASS(I) + SUM 
20 CONTINUE 
XPB(J) = SUM 
10 CONTINUE 


RETURN 

END 
KAKKKKKKARRKAKKAKKKKAKKEKAAKARKRRKKKKAAK AAA KKARKKKKKK KK 
- MODATB: READ AND TAYLOR MODAL DATA OF BODY B x 


KAKKKKKKKKARRKARAKAKAKKAAKEKARRARRARKRKAKAAKKKAKAAKAARKKRK KKK K 
SUBROUTINE MODATB (M,NPT, Xl) x2) x3 ,PISx ieee 


= PS1,PS2,PS3,CMASS , FREQ, OMEGA) 
C 
IMPLICIT REAL*8(A-H,0-Z) 
DIMENSION FREQ(M),OMEGA(M),CMASS (NPT) 
DIMENSION X1(NPT),X2(NPT) ,X3(NPT) 
DIMENSION FIBX(NPT,M),FIBY(NPT,M),FIBZ(NPT,M) 
DIMENSION PS1(M),PS2(M),PS3(M) 
DIMENSION PSI1(100,30) ,PSI2(100,30) ,PSI3(100,30) 
DIMENSION VB(3),VC(3),V0(3),VN(3) ,C(3,3) 
G 
REWIND 8 
REWIND 21 
C 
Pl) 4 0m en Ar een en) 
el 2.6) iu 
FRAD = PI/180.0 
€ 
READ(8,100) ICOR,IMAS,ICMOD,ICTR,NGRID,NMOD 
€ 
TPCECIR 2 LE mGommeme 
IF(ICTR GE. 3) (Gomnon 
READ(8,104) ALP 
GO To 3 
2 READ(8,106) VB(1),VB(2),VB(3), 
+ VC(1) VWel20e Ves) 
€ 
3 IF(IMAS .EQ. 0) GOTO 9 
c 


DO 5 I=1,IMAS 


68 


2) 


20 


30 


40 


50 


60 


65 
70 


READ(8,150) J1,J2,BMAS 
DO 4 J=J1,J2 
CMASS(J) = BMAS 

CONTINUE 


DO 10 I=1,NGRID 
Renmas, 200) x1(1) x2¢1) x31) 


TEC ICoReaeO. ©) GO 10-30 


DO 20 I=1,NGRID 

eae oles 

THTA = X2(I) * FRAD 
X1(I) = R * COS(THTA) 
K2(I) = R * SIN(THTA) 


DO 50 I=1,NMoD 

READ(8,300) INOD 

READ(8,400) FREQ(IMOD) 
OMEGA(IMOD) = TPI * FREQ(IMOD) 


DO 40 J=1,NGRID 
READ(8,500) FIBX(J,IMOD) ,FIBY(J,IMOD),FIBZ(J,INOD), 
PSI1(J,IMOD) ,PSI2(J,IMOD) , PSI3(J,IMOD) 
CONTINUE 
PS1(I) = PSI1(NPT,I) 


FS2(I) = FSI2(NPT,I) 
PS3(I) = PSI3(NPT,I) 
CONTINUE 


Paerein «LE. 1) GO TO 80 
Peter. .Gk. 3) GO TO 75 


Z2-DIMENSIONAL COORDNATE TRANSFORM 


ALP = ALP * FRAD 


CP = COS(ALP) 

SP = SIN(ALP) 

DO 60 I=1,NGRID 

x Sel (1) 

Zac) 

mee eux * CP + Z * SP 
K3(I) = - X * SP + Z * CP 


DO 70 I=1,NMOD 
DO 65 J=1,NGRID 
= esd 1) 
Z = FIBZ{J,I) 


FIBX(J,I) = peomer + 2.) sor 
ead) =X = SP + Z = CP 
CONTINUE 
GO TO 80 


3-DIHENSIONAL COORDINATE TRANSFORM 
69 


a 


O 


V2 


76 


el 
78 


80 


e3 
85 
837 


ae 


218) 
25 


100 
104 
106 
Ee 
is) 
160 
200 
20 
S00 
400 
410 
500 


CALL CERNG( VE VC7e> 


DO 76 I=1,NGRID 
VO(1) = X1(T) 
VO(2) = X2(I) 
vo(3) = X3(I) 

CALL TRN3(VO,VN,C) 


xT = Vine) 
X2(I) = VN(2) 
X3(I) = VN(3) 


DO! 76. 1=1en oD 
DO 77 iateeR ID 
VO(1) = FIBX(J,I) 
VOX 2) a=) Enea 
VO(3) = FIBZ(J,I) 
CALL TRN3(VO,VN,C) 
FIBX(J,I) = VN(1) 
FIBY(J,I) = VN(2) 
FIBZ(J,1) = VN(3) 

CONTINUE 


WRITE CUTEU ELLE 


WRITE(21,110) NGRID,NMOD 

IF(IMAS .EQ. 0) GO TO 85 

DO 83 I=1,NGRID 

WRITE(21,160) I, CMASS(I) 

DO 87 I=1,NGRID 

WRITE(21,210) 2, KICP)ex2CI)e aim 


DO 95 I=1,NMOD 
WRITE(21,410) I, FREQ(I), OMEGA(I) 
DO 90 J=1,NGRID 
WRITE(21,550) J, FIBX(J,1I), FIBY(J,1), FIBZ(J,1) 
IF(ICMOD .EQ. 0) GO To 95 
DO 93 J=1,NGRID 
WRITE(21,550) J, PSI1(J,1I), PSI2(J,1), PS13(J,1) 
CONTINUE 


FORMAT 


FORMAT (10,615) 

FORMAT (F10.4) 

FORMAT (10X,3F10.4/10X,3F10.4) 

FORMAT(/'NUMBER OF GRIDS =',I5/'NUMBER OF MODES =',I5) 
FORMAT(10X,21I5,D15.4) 

FORMAT(I4,D15.6) 

FORMAT (20X,3F10.3) 

FORMAT (2X,14,3X,3D15.6) 

FORMAT (8X,I4) 

FORMAT (8X,D14.6) 

FORMAT(/'MODE =',I5/'FREQUENCY =',F10.3,2X,'HZ'/'OMEGA =',D15.6) 
FORMAT (2X,6D13.6) 
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poO)  PORHMAI(3xX,14,3X,3D15.6) 
C 
RETURN 
END 


KARKKRKKRKRRKARKK KK AK KAKAAKKKKKKAKKKKKAKKKKAKKKKKAAKAKKKAKKKK 


x MODATC: READ AND TAYLOR MODAL DATA OF BODY C  * 
KEAEKEKRKARRKKKKRKARKK EK AKRKAKKKAKKKARKKKAAKAAKKARKKKKKKKK 
SUBROUTINE MODAIC (M,NPT,X1,X2,X3,FIBX,FIBY,FIBZ, 

+ CMASS , FREQ, OMEGA) 


IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION FREQ(I{) ,OMEGA(M) , CMASS(NPT) 
DIMENSION K1(NPI),X2(NET),X3(NPT) 

DIMENSION FIBX(NPT,M),FIBY(NPT,M) ,FIBZ(NPT,M) 
DIMENSION VB(3),VC(3),VO(3),VN(3),C(3,3) 


REWLND 3 


PI = 4.0 * ATAN(1.D0) 
eee 20) % Pi 
FRAD = PI/180.0 


Reo .OC) ICOR, IMAS,ICHOD,ICTR,NGRID,NMOD 


[eEGrRete wi) GO TO 3 
MaolernR Ge. 3) Go TO 2 
READ(9,104) ALP 
Go 10.3 
2 READ(9,106) VB(1),VB(2),VB(3), 
: VOC 1) VC( Zev ers) 


3 IF(IMAS .EQ. 0) GOTO 9 


DO 5 I=1,IMAS 
READ(9,150) J1,J2,BMAS 
DO 4 J=J1,J2 
4 CMASS(J) = BMAS 
5 CONTINUE 


9 DO 10 I=1,NGRID 
10 Rese. 200) X1(1),K2(1) ,X3¢1) 


PereoR @EO. 0) GO TO 30 


DO 20 I=1,NGRID 

mS I) 

Tegee= 201) ~ FRAD 

ee = he ees THTA) 
20 X2(I) = R * SIN(THTA) 


30 DO 50 I=1,NMOD 
READ(9,300) IMOD 
READ(9,400) FREQ(IMOD) 
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OQ 


(2) 


40 
50 


60 


65 
19) 


qT 


76 


77 
fee 


OMEGA(IMOD) = TPI * FREQ(IMOD) 


DO 40 J=1,NGRID 
READ(9,500) FIBX(J,IMOD) ,FIBY(J,IMOD) , FIBZ(J,IMOD) 
CONTINUE 

CONTINUE 


IF (ICTR he eGom Ome 0 
IF(ICIR -Ge. cpeconton7> 


2-DIMENSIONAL COORDNATE TRANSFORM 


ALP = ALP * FRAD 

CP = COS(ALP) 

SP = SIN(ALP) 

DO 60 I=1,NGRID 

x Sao) 

@ = X3(I) 

ON) SO ates 45 A oe ie 
X3(D) == K * SPe+ Z * GP 


DO 70 I[=1,NNOD 
DO 65 J=1,NGRID 
X = FIBX(J,I) 
Ge FUBZ 31) 


PERC n) =) <aGreec ence 
FIBZ(J,I) =- X * SP + Z * CP 
CONTINUE 
GO TO 80 


3-DIMENSIONAL COORDINATE TRANSFORM 
CALEMEERNS (VB VC Ae) 


DO 76 I=1,NGRID 


VOC ie = cna) 
VOX Zee ai) 
VO(3) = X3(I) 
CALL TRN3(VO,VN,C) 
Xi(I) = VN(1) 
X2(I) = VN(2) 
x3(I) = VN(3) 


DO 78 I=1,NMOD 
DO 77 J=1,NGRID 
VO(1) = FIBX(J,I) 
VO(2)) =) ETEoen) 
NOs = a2) 10) 
CALL TRN3(VO,VN,C) 
FIBX(J,1) = VN(1) 
FIBY(J,I) = VN(2) 
FIBZ(J,I) = VN(3) 

CONTINUE 


toi ot 


tou ou 


Vie OUR E Uiierh Bie 
Te 


80 


es 
85 
87 


100 
104 
106 
EG, 
150 
160 
200 
210 
300 
400 
410 
500 
950 


WRITE(21,110) NGRID,NMOD 

IF(IMAS .£0. 0) GO TO 85 

Domes i NGRID 

WRITE(21,160) I, CMASS(TI) 

DO 87 I=1,NGRID 

WemmE2l 210) 1) X1(1), x2e1) x3 (1) 


DO 95 I=1,NMOD 
WRITE(21,410) I, FREQ(I), OMEGA(I) 

DO 90 J=1,NGRID 

WRITE(21,550) J, FIBX(J,I), FIBY(J,I), FIBZ(J,I) 
CONTINUE 


FORMAT 


FORMAT (10X,6I5) 

FORMAT (F10.4) 

FORMAT (10%,3F10.4/10X,3F10.4) 

FORMAT(/'NUMBER OF GRIDS =',I5/'NUMBER OF MODES =',I5) 
FORMAT(10X,2I5,D15.4) 

FORMAT (I4,D15.6) 

FORMAT (20X,F10.2,F10.4,F10.7) 

FORMAT(3%,14,3%,3D15.6) 

FORMAT (8X,14) 

FORMAT (8X,D14.6) 

FORMAT(/'NMODE ='!,15/'FREQUENCY =',F10.3,2X,'HZ'/'OMEGA =',D15.6) 
FORMAT (24X,3D15.6) 

FORMAT (3X,14,3X,3D15.6) 


RETURN 
END 


REAR KR KAKA KRKREARKRRE RAR KR RKAAKKKRKKRRKKRKRREKRKRRKRKRERAEKRKE 


* 


10 


20 


SUEROUDINE CIRN3 * 


SUBROUTINE CTRN3 (VB,VC,C) 
IMPLICIT REAL*8(A-H,0-Z) 
[meee 1,12,13,1P2,1P3,J3;K 
DIMENSION vB(3),VC(3),C(3,3) 
DIMENSION IP2(3),IP3(3) 
DIMENSION ¥N(3),YN(3),ZN(3) 
Deer? /2,3,1/ 

DATA IP3 /3,1,2/ 


BEe= 0), 

Domlomi=i.3 

BL = BL + VB(I) * VB(I) 
BL = DSQRT(BL) 

0) 2) Hes 

ZN(I) = VB(I)/BL 

(1,3) = ZN(I) 


Domaoul=1, 3 
i= iP 21) 
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30 


40 


50 


60 


I3 = 1P3() 

YN(I) = ZN(12) * VC(I2y = 2N(ia aero 
i = Oy 

ae 20) Weal © 

YL = YL + YN(I) * YN(I) 

YL = DSORT(YL) 

DO 50 T=1.2 

YN(I) = YN(I)/YL 

Cl 2s = oem) 


DO 60 I=1,3 
I2 = IP2(1I) 
3 = es) 


RN(T) >= YN(1I2) * ZN(1I3) = YN(13)) ez Ger 
C(I,1) = XN(I) 

RETURN 

END 


KRARKRKRKKKKKKRRKKRRKERKEKRRRRKKRKAKKKRRARAKARRKRRRKRAARAKKKKRKK 


x 


SUBROUTINE TRN3 us 


KRARKKKAAKAKKKKRRKARKKKRKKKRKRKRKKRKKRAKAKRKRKKRRARREREKRKKKRKAKK 


18, 


SUBROUTINE TRN3 (VO,VN,C) 
IMPLICIT REAL*8(A-H,0-Z) 
INTEGER I,J,K 

DIMENSION VO(3),VN(3),C(3,3) 


DOVIiGel=l.2 
VN(I) = 0. 
DO 10 J=1,3 
VN(I) = VN(I) + C(I,J) * VO(J) 
REGURN 
END 
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OC ae. 
PARAM 


Do Oo OSes 


a 


ARRAY 
ARRAY 
ARRAY 
ARRAY 


CALCULATE A MATRIX ( UDOT COEFFICIENTS) 

IER, IPVT,NPT,M,N,NPTC,MC 

HET tirTe =63 M=4 lic=4 

DIMENSION aA(9,9 

DINENSION U1D21(10,10),U1D22(10,10) ,U1D32(10,10) 

DINENSION UJD21(10,10),UKD21(10,10) ,BU222(10,10) ,BU223(10,10) 
DIMENSION B22UUJ(10,10),B22UUK(10,10),B32UUJ(10,10),U1D31(10,10) 
DINENSION B3U213(10,10),B32UUK(10,10) 

DIMENSION FIBX(17,10),FIBY(17,10),F1BZ(17,10),PSI1(10) 
DIMENSION FICX(63,10),FICY(63,10),FICZ(63,10),PS12(10) 
DIMENSION PSI3(10),X1(63),%2(63) ,%3(63) 

IPVT(20) 

U1D12(10) ,BU221(10) ,U1D13(10) ,UJD11(10),UKD11(10) 

B12UUJ(10) ,B12UUK(10) ,B3U211(10) 
STRIUB(10),STRNC(10),OMGAB(10) ,OMGAC(10) 


Pie Genet) 1s. L 


FARAH 


Clee e2-10,,C3=24 65, DELT=.010 


INITIAL 


* 


O) 


eI) 
gs 
1Q2 


READ IN VALUES OF CONSTANTS 

READ (25,107)IDUN 

BEAD) (25,102) U1DI11 

READ (25,107)IDUI 

DO 1 J= 1,10 

READ(25,103)L,U1D12(J) ,UJD11(J),B12UUJ(J) ,BU221(J),OMGAB(J) 
CONTINUE 

READ (25,107)IDUM 

DO 2 K = 1,10 
READ(25,103)L,B12UUK(K) ,U1D13(K) ,UKD11(K) ,B3U211(K) , ONGAC(K) 
CONTINUE 

PEAD (25,107)IDUM 

Dome = 1,10 

DOutes =41 10 

Emi comioa ts) Ujpeil ty uinei(s 1) bU222(J,1),B22UUI(J,1) 
CONTINUE 

CONTINUE 

READ (25,107) IDUM 

Nome a1, 10 

10 Get 91,10 
PEAD(25,104)L,L,U1D32(I,K) ,B3U213(1,K) ,B32UUK(I,K) 

CONTINUE 

CONTINUE 

READ (25,107)IDUM 

DO 7 K =1,10 

BO fie = ab als 
READ(25,106)L,L,U1D22(J,K),UKD21(J,K) ,BU223(J,K) ,B22UUK(J,K) 
CONTINUE 

CONTINUE 

READ (25,107)IDUM 

De CR ar ie 

Do SS) = ah ale 

READ (25,105) L,L,U1D31(J,K),B32UUJ(J,K) 

CONTINUE 

CONTINUE 

FORNAT (D15.6) 


(ES 


iS 
104 
LOS 
Le 
107 


a5 
od 


a7 
96 


eke: 
32 


900 


elit 
DERIVA 
NOSORT 


x 


10 


20 


FORMAT (14,5D13.6) 
FORMAT(214,3D15.6) 
FORMAT (214,2D15.6) 
FORMAT (214,4D15.6) 
FORMAT (15) 
DO 94 J = 1,10 
DO 95 I =1,NPT 
READ(25,900)FIBX(I,J),FIBY(1I,J),FIBZ(I,J) 
CONTINUE 
CONTINUE 
DO 96 J = 1,10 
DO 97 I =—tenene 
READ (25,900) FICX(1T, J), FICY(1, a) pe neaaneuD 
CONTINUE 
CONTINUE 
DO Sa). —eietO 
READ (25,900) PSI1(J),PSI2(J),PSI3(J) 
CONTINUE 
pO 92 I =1,NPTC 
READ (25,900) X1(I),X2(I),X3(I) 
CONTINUE 
FORMAT (3D20.6) 
CALL STRNEG (M,OMGAB,STRNB) 
CALL STRNEG (MC,OMGAC,STRNC) 
DO 91 J =1,M 
WRITE (21,900) PSI1(J),PSI2(J),PSI3(J) 
CONTINUE 
TIVE 


TA = 10.-10.*STEP(6.85) 

‘ele 2D. 

TAOL = 80.*TIME**2 

TAO2 = 320.*(1.-STEP(S.)) 

TA = SWITCH(TIME .LE. T1,TAOL,TAOZ) 
TA = 5. 

TA + 10.-10.*STEP(6.85) 


A110B = 0. 
B12UJ = 0. 
DO 40 3o= 99m 


AlIOB = O(U+1) * UlDE2 Ge stoR 
BI2ZUJ = BIZUUI(I) 4 Ud) + Be 
DO 10h =3 7 
A(1+J,1+I1) 
Cae Sl 
CONTINUE 
Aion aor 
BJU2ZQB = O. 
BJ 2ZUUJ a 
DO” 20 ie = he 
A1J10B =U1D21(J,1)*Q(1+1I) + A1J10B 
BJU2Z0B= BU222(J,1)*Q(I+1) + BJU2QB 
BJZUUJ=BZ22UUJ(d 2) “UC 1s 1 eee Zee 
CONTINUE 
Mule = O- 


Ui ete cleale: 
A le tla) 


i 


} 


76 


30 


40 


op, 


60 


200 


300 


BIUZOC = 0. 
BIZUUK —eon 
DO 30 K = 1,MC 
ALJIOCG= UlD2Z2(3,K) *O(1+hFK) + AlJ1QC 
A(1+H+K,1+J)=UKD21(J,K) 
Nee ake ee niniC t a) 
BJU2QC= BU223(J,K)*Q(14+H+K) + BJU2QC 
BJ2UUK = B22UUK(J,K) *U(1+M+K) 
CONTINUE 
A(1+J,1) = UJD11(J) + A1J10B + A1J19C 
A(1,1+J) = A(1tJ,1) 
B(1+J)=U(1)*U(1)*(BU221 (J )+BJU2ZQB+BJU2ZQC)... 
+2*U(1)*(BJ2UUJ + BJ2UUK)- STRNB(J)*Q(1+J) 


CONTINUE 

Poe =. 0. 

B12UK = 0. 

DO 60 K = 1,NC 

B12UK = Bl2UUK(K)*U(1+M+K) + B12UK 
FinghOe <= 0). 

ALK1OB = O. 

BKU20C = 0. 

BA2UUK = 0. 


BOVSsO0 f£ = 1 Hc 
ALK1QC = QO(14+M+I) * U1D32(1I,K) + A1K10C 
BKU2Q0C = B3U213(1,K)*Q(1+MN+I) + BRKU2QC 
BK2UUK =B32UUK(I,K)*U(1+M+I) + BK2UUK 
ACI+H+I,1+H+K) = 0. 


CONTINUE 
BKU20B = 0. 
BK2UUJ = 0 


DOe5> ed. = 15H 
A1K10B = Q(1+J)* U1D31(J,K) + A1K10B 
BKU2ZOB = BU223(J,K) *O(1+J) + BKU2QB 
RK2UUJ =B32UUJ(J,K)*U(1+J) + BK2UUJ 
CONTINUE 
Al10C = Q(1+M+K) * U1D13(K) + A119C 
Preeti, t+MtK) = 1. 
A(1+M+K,1) = UKD11(K) + A1K1QB + A1K10C 
eek os AC ItMGK, 1) 
B(1+M+K) = U(1)*U(1)*(B3U211(K)+BKU2QB+BKU2QC) + 
2*U(1)*(BK2UUJ+BK2UUK) -STRNC(K)*Q(1+M+K) 
CONTINUE 
A(1,1)= U1D11 + A110B + A110C 
B(1) = 2*U(1)*( B12UJ + B12UK )+ TA 
WRITE(22,200)((A(I,J),J=1,1+M+MC) , I=1,1+M+MC) 
FORMAT(5E14.6) 


CALL DGEFA(A,9,9,IPVT, IER) 
WRITE(8,300) IER 

FORMAT (120) 

IF(IER.NE. 0) GO TO 100 

Gnu DGESE (2 9.9, 1PVT,B 0) 
WRITE(4,200) (B(I),I=1,1+M+tMC) 
U =INTGRL(O.,B,9) 


Th 


Q = INTGRL(0.,U,9) 
“CALCULATE VARIOUS DISPLACEMENTS AT THREE LOCATIONS AND ROTATIONS AT TIP 
FHDISP = sor 
YHDISP = QO. 
ZHDLSP = CG, 
AHROT = QO. 
YHERO? =—0: 
ZHROT = Q. 
DO 703i = 1 
XHDISP = KHDISP + Q(I+1)*FIBX(NPT,I) 
YHDISP = YHDISP + O(I+1)*FIBY (NPE 
ZHDISP = ZUDISP + OC I+])-FIBZCicn ae 
AHROL ="AHROn + OCl tier os ee 
YHROT = YHROT + Q(I+1)*PSI2(I) 
ZHROT = ZHROT + Q(I+1)*PSI3(I) 
70 CONTINUE 
AUREL 0% 
YLREL = Q. 
4URE DL = 06 
SOR OF 
Y4REL = OQ. 
ZERE Ee =O. 
DO SU =a MG 
XLREL = XLREL + Q(I+1+M)*FICX(1,I) 
YLREL = YLREL + Q(I+1+M)*FICY(1,1) 
ZIREL = ZIREL + QO ( Pie) 46 CZ ieee: 
K4REL = X4REL + Q(I+t1+M)*FICK(4,1I) 
Y4REL = Y4REL + Q(I+14+M)*FICY(4,1) 
Z4REL = Z4REL + Q(I+1+M)*FICZ(4,1) 
30 CONT BCE 
S1IDISP = XHDISP + X1REL + (YHROT*X3(1) - ZHROT*X%2(1)) 
YIDISP = YHDISP + YLREL + (ZHROT*X1(1) - XHROT*X3(1)) 
ZIDISP = ZHDISP + Z1REL + (XHROT*X2(1) YHROT*X1(1)) 
X4DISP = XHDISP + X4REL + (YHROT*X3(4) - ZHROT*X2(4)) 
Y4DISP = YHDISP + Y4REL + (ZHROT*X1(4) XHROT*K3(4) ) 
Z4DISP = ZHDISP + Z4REL + (XHROT*X2(4) - YHROT*X1(4)) 
RETURN 
100 WRITE(G, 101) Tes LER 
OM FORHAT('O IER =',17) 


CALL ENDJOB 
PRINT TA,Q(1-5),U(1-5) 
CONTRE FINTIM=TO  PEEERT aul 
SAVE OL > parcels) HOLS) pete) ,Q9(4) PONS) ,9(6) DONT RCS) oD) mea a! 
U(1),XHDISP, YEDISP ,ZHDISP  AHROn, YHROT=7hRGt up Loree. 
YIDISP, 2ZUDIsr  A40 LS eye bse eee 


XRAPH (G1,NI=7,DE=TEK618) TIME(NI=5,LE=10.,UN= 

* 'SEC'), TA(UN='LB-IN') 

*RAPH (G2,NI=7,DE=TEK618) TINE(NI=5,LE=10.,UN='SEC' ) 
x ,0(1)(UN='RAD') 

*RAPH (G3,NI=7,DE=TEKG18) TIME(NI=5,LE=10.,UN='SEC' ) 
x ,U(1) (UN=!RAD/SEC'! ) 

*RAPH (G4,NI=7,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC' ) 
0(2)(UN='IN') 

*PAPH (GS,NI=7,DE=TEK618) TINE(NI=5,LE=10.,UN='SEC') 
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ms 
RE 


~ 


*RAPH 


~ 


*RAPH 
* 


Saal 
*& 


“PAPH 
* 


*RAPH 
~ 


*RAPH 


x 


GRAPH 


*PAPH 


ras 


GRAFH 


*RAPH 


Bas 


*RAPH 
* 


*RAPH 
* 


GRAPH 


, 0(3)(UN='IN') 

(G6 ,NI=7,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC' ) 

(0(4) (UN="III') 

(G7 ,NI=7 ,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC') 

,O(5)(UN='IN') 

(G7A,NI=7 ,DE=TEK618) TIME (NI=5 ,LE=10.,UN='SEC' ) 
(0(6)(UN='IN') 


(G7B , NI=7 ,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC') 
OOS ie 
(G7C, NI=7,DE=TEK613) TIME (N1=5,LH=10.,UN='SEC’ ) 
,O(8) (UN="IN") 
(G7D,NI=7,DE=TEK618) TIME(NI=5 ,LE=10.,UN='SEC' ) 


, Q(9)(UN='IN') 

(G8,NI=7,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC') 
PIG r UN Mili hl=5,LO=—10. St-2. } 

(G9,NI=7 ,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC') 
Dio? ( UN wi LI=4,10=-2, ,SC=1a) 

(GA ,NI=7,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC') 
/ZHDISP(UN='IN! ,LI=4) 

(GR8,'NI=7,DE=TEK618) TINE(NI=5,LE=10.,UN='SEC') 
/XHROT(UN='RAD' ,LI=4,LO=~.020,SC=.005) 
(GR9,NI=7,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC') 
/YHROT(UN=!RAD') 

(GRA ,NI=7,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC') 
/ZHROT(UN='RAD' ) 

(GB,NI=7,DE=TEK618) TIME(NI=5,LE=10.,U!='SEC' ) 
Pam Cui" Mi ml=5 Le=-20. ,SC=5. ) 

(Gemiii=7 De=PHkel6) TIME(NI=5 ,LE=10. ,UN="SEC') 
PalpoSr U='IN" Li]=4,1L0=-3.,SC=2. ) 


*GRAPH (GD,NI=7,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC' ) 


¢ 
-_— 
as 


Pei DLSPC(UNS'IN"') 


*GRAPH (GE,NI=7,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC' ) 


oe 


GRAFH 


Perse CUN='IN',Li= 5,L0=-8. ,SC=2.) 
(GF ,NI=7,DE=TEK618) TIME(NI=5,LE=10.,UN='SEC' ) 
/Y4DISP(UN='IN' ,LI=4,L0=-1.2,SC=1.) 


*GRAPH (GG,NI=7,DE=TEK618) TINE(NI=5,LE=10.,UN='SEC' ) 


LABEL 
LABEL 
L&BEL 
EASE |, 
LABEL 
LABEL 
LABEL 
LABEL 
LABEL 
Peek 
LABEL 
LABEL 
Eyete 
LABEL 
Easel 
LABEL 
LABEL 


,Z4DISP(UN='IN') 
(G1) APPLIED TORQUE 
(G2) ANGULAR DISPLACEMENT 
(G3) ANGULAR VELOCITY 
(G4) GENERALIZED DISPLACEMENT Q1 
(G5) GENERALIZED DISPLACEMENT 92 
(G6) GENERALIZED DISPLACEMENT Q3 
(G7) GENERALIZED DISPLACEMENT 04 
(G7A) GENERALIZED DISPLACENENT 95 
(G7B) GENERALIZED DISPLACEMENT Q6 
(G7C) GENERALIZED DISPLACEMENT 97 
(G7D) GENERALIZED DISPLACEMENT 98 
(G3) X DEFLECTION AT BOOM TIP 
(G9) Y DEFLECTION AT BOOM TIP 
(GA) Z DEFLECTION AT BOOM TIP 
(GR8) X ROTATION AT BOON TIP 
(GER9) Y ROTATION AT BOOM TIP 
(GRA) Z ROTATION AT BOOM TIP 
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ee ¢ ¢@ 


LABEL (GB) 
LABEL (GC) 
LABEL, (ED) 
LABEL (GE) 
LABEL (GF) 
LABEL (GG) 
END 

STOP 
FORTRAN 
KAKKKKKRKKKKR RR RK KR KRKKKRK RAK KEK KKK KR KKK KR EK KARR KKK KK KKK KKK KKKKKKKKRKKKKKKK 


x STRAIN ENERGY S 
KEE RAKA RRR KIRK K RK KK RRR KR KKK RK RRR ERR KEK RRR KR KKK RK KKK KKK KKK KK KKK AKAAKH 


DEFLECTION At Pish seein 
DEFLECTION TAL Dis7 soll 
DEFLECT FONeAd Dish eOot 
DEFLECT TON. Ss) DISh FOLia 
DEFLECTION At Disheeeina 
DEFLECTIONPAL Diseer one 


Nir oS NOK XX 
PoP RR Re 


SUBROUTINE STRNEG (M,OMEGA,STRAIN) 
REAL*8 OMEGA(10) ,STRAIN(10) 
DOM Oe ane 
STRAIN(I) = OMEGA(I)**2 
10 CONTINUE 
RETURN 
END 
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Figure 2.1 


NROSS Baseline Configuration 
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LFMR Model 
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Figure A.1 


Undeformed Reflector 
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Figure A.6 


Undeformed Boom 
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Figicees a7 


First Mode of Boom 
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Figure A.8 


Second Mode of Boom 
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Figure A.9 


Third Mode of Boom 
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Fourth Mode of Boom 
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Figure 4.1 
Finite Element Model of LFMR System 
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Applied Torque History During the Spin-up Procedure 
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Figure B.1 


Angular Velocity of Stiffer Reflector at -155° 
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Figure B.2 


Angular Velocity of More Flexible Reflector at -155° 
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Fourth Boom Generalized Coordinate of Stiffer Reflector at -155° 
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Fourth Boom Generalized Coordinate of 
More Flexible Reflector at -155° 
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Third Reflector Generalized Coordinate of 
Stiffer Reflector at -155° 


100 


TIME (SEC) 


10 





e e e e e e 
a w w 
o = a Nw r 6 4 


(NT 6-06) LO 1IN3W3901dS10 0321 WYSN39 


Figure B.6 


Third Reflector Generalized Coordinate of 
More Flexible Reflector at -155° 
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Figure B.7 


Vertical Deflection at Boom Tip of Stiffer Reflector at -155° 
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Figure B.8 


Vertical Deflection at Boom Tip of 
More Flexible Reflector at -155° 
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Figure B.9 


Rotation of Boom Tip About Vertical Axis 
of Stiffer Reflector at -155° 
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Figure B.10 


Rotation of Boom Tip About Vertical Axis 
of More Flexible Reflector at -155 
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Figure C.1 


Angular Displacement of Stiffer Reflector at -135° 
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Angular Displacement of Stiffer Reflector at -145'° 
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Figure C.3 


Angular Displacement of Stiffer Reflector at -155° 
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Figure C.4 


Horizontal Deflection of Boom Tip of Stiffer Reflector at -135° 
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Horizontal Deflection of Boom Tip of Stiffer Reflector at -145' 
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Figure C.6 


Horizontal Deflection of Boom Tip of Stiffer Reflector at -155° 
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Rotation of Boom Tip About Horizontal Axis 
of Stiffer Reflector at -135° 
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Figure C.3 


Rotation of Boom Tip About Horizontal Axis 
of Stiffer Reflector at -145° 
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Figure C.9 


Rotation of Boom Tip About Horizontal Axis 
of Stiffer Reflector at -155° 
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Horizontal Deflection of Dish Point 1 of Stiffer Reflector at -135'° 
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Horizontal Deflection of Dish Point 1 of Stiffer Reflector at -145° 
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Horizontal Deflection of Dish Point 1 of Stiffer Reflector at -135° 
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Figure D.1 


First Boom Generalized Coordinate of Stiffer Reflector at -155° 
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First Boom Generalized Coordinate of 
Stiffer Reflector Tilted 5° Out of Plane 
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Figure D.3 


Fourth Reflector Generalized Coordinate 
of Stiffer Reflector at -155° 
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Figure D.4 


Fourth Reflector Generalized Coordinate 
of Stiffer Reflector Tilted 5° Out of Plane 
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Out-of-Plane Deflection at Boom Tip of Stiffer Reflector at -155° 
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Figure D.6 


Out-of-Plane Deflection at Boom Tip of 
Stiffer Reflector Tilted 5° Out of Plane 
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Out-of-Plane Deflection at Dish Point 4 
of Stiffer Reflector at -155° 
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